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A Fundamental Comparison of Incomplete 
Charge Transfer in Charge Transfer 


Devices 


By C. N. BERGLUND*® and K. K. THORNBER 
(Manuscript received September 22, 1972) 


UJsing a small-signal analysis, we present a general comparison of the 
more tnportant contributions to the ancomplete transfer of charge in the 
following charge transfer devices: (2) the polyphase charge-coupled device 
and the [GIF ET bucket-brigade shift register, where individual charge 
transfers are single-step processes, and (221) the two-phase charge-coupled 
device, the conductively connected charge-coupled device, the tetrode bucket 
brigade, and the stepped-oxide bucket brigade, where individual charge 
transfers are two-step processes. A recently proposed lumped-charge-model 
approximation 1s made in order to estimate the teme dependence of the 
transferred charge including both drift and diffusion. In this calculation 
we also include the effects associated with the injection of charge from a 
diffused source into the IGFET channel which modify the current-voltage 
behavior at low currents. Using this calculation of the tume dependence of 
the transferred charge, the various contributions to incomplete transfer, 
encluding those due to trapping in the interface states, are derived and 
compared for each of the charye transfer devices of interest. The results 
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show that highest transfer efficiencies at low frequencies will usually be 
obtained from suitably designed shift resisters utilizing a two-step transfer 
process. This ts because a two-step process tends to reduce those contribu- 
trons to incomplete transfer which are due to signal-charge-modulation of 
the conductances and capacitances governing the charge transfer current. 
For example, the coefficient of incomplete charge transfer, a, for an indi- 
vidual transfer in a polyphase CCD is found to be 0.6, 0.85, and 1.4-10-% 
at 2, 4,and 8 MHz, respectively. By contrast, for the same physical device 
dimensions, a in a stepped-ovide, two-phase CCD 1s reduced to 0.38, 0.6, 
and 0.8-10~ at the same frequencies, nearly a factor of two reduction. 
The intrinsic contributions to incomplete charge transfer, which dominate 
the high-frequency operation, are found to be relatively insensitive to the 
type of device considered for the dimensions given. An a of 0.8-10~ at 
100 MHz ts typical. Among the various two-step process devices, interface- 
state effects are found to be comparable, 0.3-10~% being a typical value for 
this frequency-independent contribution. At low frequencies, a frequency- 
independent contribution to incomplete transfer in bucket-brigade transfer 
processes 1s found. The size of this contribution 1s about 0.15-107%. These 
results further illustrate the utility of analyzing incomplete transfer effects 
mm terms of single-device, small-signal characteristics. 


I. INTRODUCTION 


Since the introduction of the charge-coupled device! (CCD) and the 
bucket-brigade shift registers,?~* there have been a number of proposed 
modifications to the basic structures in order to achieve specific per- 
formance goals.*- It is the purpose of this paper to analyze the 
mechanisms of charge transfer in these devices and to compare the 
characteristics of the incomplete charge transfer in each. In this way 
one can evaluate the relative merits and limitations of these modified 
charge transfer devices (CTD’s). 

In a recent article, the authors!! presented a general analysis of in- 
complete charge transfer in CTD’s. Using a lumped model”?:® to 
characterize the dynamics of the charge transfer, it was possible to 
calculate a, the small-signal coefficient of incomplete transfer, in terms 
of single-device, small-signal characteristics. It was found that three 
contributions to incomplete transfer are common to all charge transfer 
shift registers: an intrinsic transfer rate contribution, an output con- 
ductance or feedback contribution, and a storage-capacitance modu- 
lation contribution. A significant feature of that calculation was that 
each of these contributions could be expressed analytically in terms of 
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the charge transferred as a function of time. Analytic results are, of 
course, very convenient for comparative purposes. In addition, numeri- 
cal values for a could be obtained, the accuracy of which depends only 
on the accuracy to which the transferred charge is known, and not 
on the much lower accuracy obtainable from the difference between the 
transferred charge for slightly different signal sizes. Although originally 
developed to treat the standard CCD and bucket-brigade?* devices, 
it is possible by a straightforward generalization to analyze the more 
complicated device structures considered here with little additional 
effort. 

In this paper we derive in a uniform manner the incomplete transfer 
properties! of the simple bucket brigade, the three-phase CCD, various 
two-phase CCD’s, the conductively connected CCD, and the tetrode 
and stepped-oxide bucket brigades. These calculations are then used 
as a basis for comparison of the performance to be expected from each 
CTD. It will also become clear that other CTD’s can be treated in a 
similar fashion. 

Using a lumped-charge model for CTD behavior® based on a charge- 
control concept proposed by Lee and Heller,” we first calculate the 
time-dependence of the transferred charge in a CTD for both a bucket- 
brigade cell and a CCD cell. This derivation of necessity 1s somewhat 
approximate, but it is found to be satisfactory for the purpose of calcu- 
lating the incomplete transfer properties by our method and is com- 
parable with the other approximations necessary for any such calcu- 
lation. Being simple, it emphasizes the important characteristics of the 
transfer event. From the time dependence of the charge transfcr we 
determine the small-signal transfer inefficiency" and calculate the 
clock-frequency dependence of the several contributions to incomplete 
transfer. Finally, these theoretical results are applied to several 
specific CTD’s and their important operating limitations are discussed. 


II. TIME DEPENDENCE OF TRANSFERRED CHARGE 
2.1 Preliminary Considerations 


The starting point in determining the performance characteristics 
of a charge transfer device is a calculation of the time dependence of 
the charge transfer process. However, the nonlinearity of the problem, 
the requirement that several assumptions be made concerning the 
charge transport and trapping along a semiconductor interface, and 
the dependence of the solution on clock-voltage waveform and on the 
details of device design, all combine to make a general solution valid 
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for all cases an unrealistic goal. Indeed for our purposes here a highly 
accurate solution 1s not necessary. We shall, therefore, make assump- 
tions to simplify the problem considerably. 


2.2 General Approach 


Our first assumptions are that the charge transfer in the more com- 
plicated modified structures can be treated as a series of simpler 
transfer steps, and that each charge-transfer step can be characterized 
by either a ‘“‘CCD” transfer process or a “‘bucket-brigade”’ transfer 
process. The major difference between the two is that in the ‘‘bucket- 
brigade” transfer process the injection of carriers from a diffused 
region into a surface channel is taken into account. In the “CCD” 
transfer process, by contrast, the charge carriers are initially in the 
surface “‘channel” (storage well) and no injection over a barrier is 
necessary. We shall solve for the charge transfer as a function of time 
for both types of transfer processes. Such solutions will be sufficient to 
discuss CT D’s in which the charge transfer from one storage region to 
the next are single-step processes. For CTD’s with more complicated 
charge transfer mechanisms involving multiple-step processes, the 
charge transfer is represented as a series of single-step processes of 
either the “CCD” or the “bucket-brigade”’ type as is appropriate. 
In this way, sufficiently accurate results will be achieved to illustrate 
most of the important features of the charge transfer and to provide 
a basis for comparison of the several CTD schemes. 

In order to carry out the above program, we must calculate the 
charge transfer in a ‘“‘CCD” and in a “‘bucket-brigade”’ transfer process. 
The general approach that will be used will be the charge-control or 
lumped-charge model. This type of model has usually been used for 
time-dependent analyses of circuit problems'*—!* as well as for analysis 
of some CTD’s such as the IGFET” and bipolar!’ bucket-brigade shift 
registers. By comparison of results from such a model with exact com- 
puter solutions,!®~*! it has been shown that the model is also applicable 
with surprising accuracy to the charge-coupled device.’*? We shall 
assume that the gradual channel approximation 1s valid for all the cases 
to be treated here. This means that the one-dimensional solution to 
Poisson’s equation will be used throughout and that fringing electric 
fields will influence only a small fraction of the active device area, and 
will, by definition, mean that our results will not be applicable to CTD’s 
such as the buried channel CCD,” except possibly in the high-frequency 
limit. For simplicity, it will also be assumed that we are dealing with 
p-channel devices with a constant carrier mobility, that depletion- 
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Fig. 1—IGFET gate region illustrating the symbols defined in the text. 


layer capacitances are constant, voltage-independent quantities, and 
that the clock voltage driving the CTD’s are ideal square waves. 

We begin by calculating the steady-state transport under given 
boundary conditions for an IGFET gate. This is similar to the con- 
ventional IGFET solution except that we include carrier diffusion in 
addition to drift. The current-voltage relation we obtain will then be 
used in calculating the “CCD” and ‘“‘bucket-brigade” transfer 
processes. 

Referring to Fig. 1 and following Sze,?* the minority carrier charge 
density Qp(y) is given by 


Qrly) = CiV (y) (1) 


where C; is the capacitance per unit area of the oxide C,x in parallel 
with the silicon space charge Cgc, and V(y) is the silicon surface 
potential with reference to the surface potential value in the absence 
of minority carriers. A lateral electric field can only exist in the gradual 
channel approximation because of a gradient in Q,, and will be 
represented as — (1/C;) (dQ,/dy). Given that the current flowing at 
any value of y must be a constant J, and that current density J is 
given by 
dp 
J = qupk — qd — (2) 
dy 
where » and D are the carrier mobility and diffusion constants re- 
spectively (related by the Einstein relationship), p is the hole density, 
and E is the electric field, one can write for the current flowing from 
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the source to the drain: 





Zu —| dQ p(y) 3) 


where Z is the channel width. Placing the boundary condition that the 
surface potential V on the “‘drain”’ side of the channel is approximately 
zero (i.e., assume the drain voltage is sufficiently negative that we are 
operating in saturation) and on the ‘‘source”’ side is assumed equal to 
a voltage Vu, one obtains 


B kT 
l=-YV, V,zt+2— (4) 
2 q 
where 6 = ZuC;/L, is the conventional IGFET gain factor, and L, 
is the effective channel length. This expression differs from the usual 
saturation IGFET current expression (which ignores diffusion) only 
by the additive term 2k7'/q in the last factor. Similarly we can calcu- 
late the total charge Q stored under the gate for a given voltage V4. 
Representing the geometrical capacitance of the gate oxide in parallel 
with the underlying silicon space-charge capacitance as C'go, 
3 kT 
Va 
2 24 
Q = 3CaoVa ap ; (5) 
Va 2 
qd 
Hence, for V4 very large or very small with respect to kT’/q, Q is linear 
in the voltage V4. Thus an effective gate capacitance Cg can be de- 
fined which is % or 4 respectively of the geometrical capacitance. For 
our purposes here, this is sufficiently close to a constant value that we 
can assume Cg is equal to #C go over the entire range of V4 of interest. 


2.3 The “CCD” Transfer Process 


In the CCD problem treated by several authors,!%—?1.24 the tran- 
sient decay of the charge stored under a single capacitor plate is cal- 
culated. With the ideal square-wave clocks assumed here, this problem 
is illustrated in Fig. 2. Given the charge-control or lumped-charged 
model assumption, the problem reduces to one of discharging a 
capacitor through a nonlinear resistor derived from eq. (4). In the 
lumped-charge model as shown above, the total charge under the plate 
is approximately linear in the voltage V4 (see Fig. 2). Hence, one can 
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Fig. 2—Transfer of charge in a standard CCD cell with symbols defined in the text. 


write from eq. (4) 


d dV kT 

WS = Og = ~1 = 2 ¥4(Vat2—) 
dt dt 2 q 

(ZL, in 6 now refers to the effective length of the storage well of the 


CCD cell). If the initial voltage at A at the beginning of transfer was 


(6) 





10° 


107! 


V/V, FOR V, = 10 VOLTS 


10-2 


1073 


NORMALIZED VOLTAGE, CURRENT 





0 1 4 9 16 25 36 49 
NORMALIZED TIME (Bt/Cg) 


Fig. 3—Normalized voltage and current plotted as a function of normalized time 
for the transfer of charge from one well to the next in a standard CCD. 
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Vand the voltage at time tis V, then eq. (6) gives the solution 


Vo qVo BkT “) 
1+ 22/1 - ex (- ‘) 
2kT qCg 

Except for the fact that V and V, are voltages at point A rather than 
some appropriately defined average voltage under the plate, this 
result is very similar though not identical to one derived by Lee and 
Heller” and is plotted in Fig. 3 assuming a 10-volt value for V.. 
Comparison of eq. (7) with the computer results of Strain and Schryer"” 
also gives excellent agreement provided that the voltage V in eq. (7) is 
scaled to reflect the average voltage as defined by Strain and Schryer.'® 
It will be shown later that in calculating the incomplete transfer 


properties of CTD’s the current J at time ¢ compared to the current 
I, at the initial part of the transfer is of importance. From eq. (6) 





V4 2kT 
V q 
— = — |—____ (8) 
Vo 2kT 
Vo + 5 a. 
q 


Oo 


and this is also plotted in Fig. 3. 


2.4 The “Bucket-Brigade’”’ Transfer Process 


The bucket-brigade transfer problem is very similar to that of the 
CCD in the sense that the current during transfer is given by eq. (4) 
with the voltage V4 being measured at point A in the channel near 
the source (see Fig. 4). The important difference which requires us to 
treat the transfer process separately is that in the bucket brigade the 
minority carriers must be injected into the channel from the diffused 
p-regions, 1.e., a large fraction of the charge is stored in the diffused 
regions and the p-island voltage is different in general from V4. 
Representing the voltage on the source p-region as V,, the storage 
capacitance associated with the p-island as C’,, and that associated 
with the channel as Cg, we can write analogous to eq. (6) 

aVa dV, B kT 
Ce +C, = -50a(Vat2—): 
dt dt 2 q 








(9) 
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Fig. 4—Transfer of charge in a standard IGFET bucket-brigade cell with symbols 
defined in the text. 


In order to solve this equation, a relation between the p-island voltage 
and V4 must be obtained. This problem has been treated somewhat 
differently by Barron.2® For our purposes here, we shall ‘simply 
assume that the minority carrier density at point A is related to the 
source p-island voltage by an expression of the form 


CiVia = Qn = Qp. expla(Vp — Va)/kT ] (10) 


where Q»y. is some constant [see eq. (1) ]. Taking derivatives with 
respect to time gives 














dV » kT \ dVa 
= (: ) (11) 
dt qV a dt 
and substituting in eq. (9) we obtain 
Vat+2 tad 
aV 4 B 2 : qd 
= —— V4 |—_—______ | (12) 
dt ZC Ce kT 
Va (1 - =| es 
Cy q 


This expression can be solved exactly for V4 as a function of time, 
and such a calculation has been carried out assuming Cg negligible 
compared to C’, and is plotted in Fig. 5. In addition to this, however, 
it is of most value to note that eq. (12) reduces to the CCD transfer 
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Fig. 5—Normalized voltage and current plotted as a function of normalized time 
for the transfer of charge from one p-region to the next in a standard IGFET bucket- 
brigade cell. V/V, from Fig. 3 is shown for comparison. 


result in the limit of Cp equal to zero, so that in that sense eq. (12) isa 
general expression for single-step charge transfer devices. Also, the 
last factor in eq. (12) for C, large compared to Cg is a very slowly 
varying function with value on the order of one or two depending on 
Va. If we assume this quantity to be unity, then the solution to eq. 
(12) is the previously* derived bucket-brigade result assuming an ideal 
square-law current-voltage relation 
Va 1 


Sit jt node Bs (13) 
Vo BV ot 
es 


2C'p 


Hence we see that for the bucket-brigade transfer process under square- 
wave clocks, the voltage V4 varies as reciprocal time after long times 
rather than exponentially as it does for the CCD transfer process. This 
is also illustrated by the exact solution shown in Fig. 5. Note from eq. 
(10) that V, can be related to V4 through 


kT Va 
Vao— Vp = ay 











(14) 
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where V,. is a constant. Thus V, follows V4 but differs by a term 
logarithmic in V4. Since V4 varies reciprocally with time at long times, 
V, ultimately has a logarithmic dependence on time as reported by 
Buss and Gosney.?® 

For comparison to the CCD transfer process, we have plotted V4 
from Fig. 3 in Fig. 5. Also of interest is the time dependence of the 
transfer current, given by eq. (8) and plotted in Fig. 5. 


III. TRANSFER EFFICIENCY IN CHARGE TRANSFER DEVICES 


Given that one has calculated the charge left behind as a function 
of transfer time as in the previous section, the problem of using these 
results to predict the performance of an n-stage shift register remains. 
Under most shift register operating conditions, a circulating charge or 
fat zero will be used, so that a large fraction of the charge left behind 
after each transfer will be independent of the signal charge. Thus, a 
definition of transfer inefficiency as the total charge remaining divided 
by the total charge to be transferred will often be of little value even 
if we had a highly accurate solution including the net amount of charge 
trapped in interface states. Further, the amount of degradation in a 
shift register will depend on the time dependence of the particular 
signal being used. Since, in addition, the charge transfer process is 
generally nonlinear, it has been found convenient to define a differential 
or small-signal parameter to characterize incomplete charge transfer 
in CTD’s; 1.e., superimpose a small-signal charge on a larger back- 
ground charge being continually transferred. It has been shown?’:?8 
that by linearizing the problem in this way, the signal degradation by 
an n-stage register can be readily characterized. For this reason we will 
focus our attention here on the small-signal incomplete transfer 
parameter a defined as the change in the charge left behind after trans- 
fer divided by the small change in signal charge, which was the cause 
of the change in charge left behind. Using derivatives 

a= sae (15) 

dQ) o 

where Q, is the total charge to be transferred. (a = a(t) since Q = Q(t).) 
In a recent paper! the authors have shown that incomplete transfer 
in CTD’s can be treated in a general way using this small-signal ap- 
proach, and that a can be separated into three components: an in- 
trinsic transfer rate term, a drain conductance or feedback term, and 
a storage-capacitance modulation term. Interface states introduce 
additional terms, but of a similar nature to these three. 
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Fig. 6—(a) General representation of a single-step transfer cell of a charge transfer 
aon (b) General representation of a two-step transfer cell of a charge transfer 
evice. 


In this previous paper! we described the charge transfer process in 
a CTD in terms of the model shown in Fig. 6a. Although applicable 
for treating two-step processes, this model is really best suited for 
single-step processes. This being the case, it is convenient to treat a 
two-step process using the model shown in Fig. 6b, which is basically 
two single-step models in series. Even more complicated charge 
transfers can be treated by including additional single-step models. 

The single-step process is illustrated in Fig. 6a using the same sym- 
bols as used previously. The charge from a storage capacitor C's is 
transferred through some nonlinear conductance to the drain capacitor 
C'p as in the simple bucket brigade or the three-phase CCD. (Vc is 
the clock voltage.) However, in the two-step process illustrated in 
Fig. 6b, the charge from the storage capacitor C's is transferred to some 
intermediate capacitor Cz then to the drain capacitor Cp during a 
single transfer period. The advantage of this is that channel-length 
modulation effects will be reduced, as will be described later. Examples 
of this kind of transfer are the two-phase CCD and the tetrode bucket 
brigade. 

In Appendixes A and B, the various contributions to the incomplete 
transfer parameter a are calculated for both single-step and two-step 
CTD?’s, and in the next section, a is evaluated theoretically for several 
proposed CTD structures in order to provide a comparison and point 
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TABLE I—List or AssuMPTIONS MADE TO SIMPLIFY 
CTD INcoMPLETE TRANSFER CALCULATIONS 


1. All devices are based on the silicon-silicon-dioxide system using n-type silicon sub- 
strates doped to a density of 10!° cm™%. 

. The silicon dioxide thickness is 1000 A. 

. All charge-storage regions, gates, and diffusions have a length of 10 microns. 

. Net interface state charge after transfer is constant at 2 x 10°° coulombs per 
square centimeter. 

. Minority carrier mobility is constant at 200 cm?/volt-second. 

. Fringing electric fields penetrate a distance equal to the one-dimensional space- 
charge width. 

. The background charge on which the signal charge is superimposed corresponds 
to a voltage at point A of 10 volts, and the drain-to-substrate voltage reaches 
a value of 10 volts at the end of transfer. 

8. Clock voltages are ideal square waves. 


~I > Or He Oo bo 


out the important incomplete transfer mechanisms. To simplify the 
calculations and provide specific examples, several assumptions have 
been made which are listed in Table I. Calculations for conditions 
other than those listed in Table I can je made using the derivations 
in the appendix. In Table IJ we summarize the nature of the transfer 
processes involved in each type of device discussed in Section IV. 


IV. INCOMPLETE TRANSFER PROPERTIES OF SEVERAL CHARGE TRANSFER 
DEVICES 


4.1 The Simple Bucket Brigade 


The simple bucket-brigade shift register is a two-phase CTD 
fabricated as illustrated in the insert in Fig. 7. In such a shift register, 
it has been shown‘! that the two dominant terms under most con- 
ditions are the drain conductance or feedback contribution, ap, and 
the intrinsic transfer rate contribution, a;. However, the analyses 


TABLE II—SUMMARIES OF THE NATURE OF TRANSFER PROCESSES 


Type of Device Mode of Charge Transfer 
Single-Step Transfer Process 
Standard CCD CCD 
Standard Bucket Brigade BB 
Two-Step Transfer Process First Transfer Second Transfer 
Two-Phase CCD CCD CCD 
C4D CCD BB 
Stepped-Oxide BB BB CCD 


Tetrode BB BB BB 


160 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1973 


1071 





INCOMPLETE TRANSFER 


0.1 0.2 0.5 1.0, 2 5 10 20 50 100 
CLOCK FREQUENCY IN MHz 


Fig. 7—Contributions to a (incomplete transfer) in the simple bucket-brigade shift 
register shown in the insert. 


have usually assumed a simple square-law IGFET current-voltage 
characteristic and have ignored the fact that the carriers must be 
injected into the channel from the diffused regions. Recently?® it has 
been pointed out that this injection requirement leads to a charge 
left behind which is logarithmic with time and an incomplete transfer 
parameter which tends to a constant value at low clock frequencies. 
Using the assumptions listed in Table I and assuming in addition that 
the drain capacitance Cp is equal to C>, Fig. 7 shows the calculated 
behavior of the two contributions to a, a; and ap, as a function of 
clock frequency. The tendency for a to saturate at low frequencies is 
apparent in addition to the previously derived linear behavior of ap 
and quadratic behavior of a; with clock frequency at the high 
frequencies. 

The other components of a, storage capacitance modulation and 
interface-state capacitance modulation, have not been shown in Fig. 7 
because their values depend on the ratio of the channel capacitance to 
the storage capacitance. Both should be relatively small compared to 
the sum of a; and ap at all frequencies, if C¢ is much smaller than C,; 
so that qualitatively the bucket-brigade shift register should be well 
approximated by considering only a; and ap. However, if C¢ and C, 
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are of the same order as will often be the case for small-area devices, 
the interface-state capacitance modulation will add a frequency- 
independent term ac,ss of value equal to approximatcly 8 X 1074 
(Cg/C,). In that case, the interface-state term may dominate a at 
low frequencies. 


4.2 The Three-Phase CCD 


As originally proposed,! the charge-coupled device was driven by a 
three-phase clocking scheme as shown in the insert in lig. 8. If it is 
assumed that the transfer time from one capacitor plate to the ad- 
jacent capacitor plate is one-third of a clock period, then the incom- 
plete transfer parameter for the CCD under the restrictions of Table I 
is calculated to vary with clock frequency as shown in Fig. 8. In 
comparison to the simple bucket brigade, the three-phase CCD has a 
similar clock frequency behavior at high frequencies and a similar 
upper limit for operation. However, at lower frequencies the two 
modulation terms ac and ac,ggs become dominant.” At very low 
frequencies, a tends to a relatively frequency-independent value due to 
the interface-state capacitance modulation term. Over most of the 
range of clock frequencies, the feedback or drain conductance term, 
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Fig. 8—Contributions to a (incomplete transfer) in the standard three-phase CCD 
shown in the insert. 
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of such importance for the bucket brigade, is nearly negligible for the 
three-phase CCD. 

In Fig. 8 the various modulation terms have been calculated assum- 
ing that the channel-length modulation can be computed using eq. (69). 
Care should be used in comparing the modulation terms in Fig. 8 with 
those in Fig. 7. Owing to the uncertainty in the estimated modulation 
of the length of the channel or storage region, the magnitudes of the 
modulation terms should be considered to be approximate. Since ap, 
ac, and ac,ss are each directly proportional to the length modulation, 
primary emphasis can be placed on their relative magnitudes and on 
the clock-frequency dependence of the various terms. 


4.3 The Two-Phase CCD 


Since the basic charge-coupled device has no inherent directionality 
like that of the bucket brigade, two-phase operation can only be 
achieved if some asymmetry is introduced into the CCD cell. One way 
to do this is to make each CCD capacitor plate consist of two regions, 
a transfer or barrier region which prevents carriers from moving in 
the wrong direction and a storage region. Ion-implanted barriers’ and 
two oxide thicknesses®? are two schemes which have been proposed, 
but two-phase operation can also be achieved by fabricating a four- 
phase CCD and placing a de bias between alternate clock lines. Figure 
9 illustrates these approaches to two-phase operation. 

It is immediately evident from Fig. 9 that a two-phase CCD operates 
using a two-step transfer process. Hence, as shown in Appendix A, 
the operational improvement over a three-phase CCD comes pri- 
marily from making C's small compared to Cp and by reducing the 
channel-length modulation. Given that the barrier length is identical 
to the storage capacitor length as assumed in this work, this means 
that best performance should be achieved either with the stepped- 
oxide device in Iig. 9, since the thicker oxide over the barrier region 
results in a smaller Cz, or with the ion-implanted device since dL./dV p 
is reduced. In calculating the incomplete transfer for the stepped- 
oxide device, we will assume that the geometrical capacitance asso- 
ciated with the barrier and the transconductance are both one-half 
the values in the storage region. 

Figure 10 shows the calculated incomplete transfer results for the 
stepped-oxide, two-phase CCD. Comparing these results to those for 
the three-phase CCD shown in Fig. 8, it is seen that the intrinsic 
transfer rate term, thus the high-frequency limitation, is about the 
same for both. The storage-capacitance modulation and drain con- 


INCOMPLETE CHARGE TRANSFER 163 


ductance terms are both reduced because of the small value of Cz, 
but the interface-state capacitance modulation term is approximately 
the same. Hence, operationally, some improvement is gained by using 
a two-phase CCD rather than a three-phase CCD. 

A word of caution concerning the results for a shown in Fig. 10 is 
in order at this point. Some recent measurements*! indicate that at 
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Fig. 9—Several examples of two-phase CCD device structures: (a) Ion-implanted 
barrier. (b) Standard four-phase run two phase. (c) Stepped-oxide. 
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ae 10—Contributions to a (incomplete transfer) for the stepped-oxide, two-phase 


low frequencies, a may be independent of interface-state density. 
However, as shown in Fig. 10, it is at low frequencies that the interface- 
state contribution to a is dominant, and this contribution is propor- 
tional to the density of interface states. The result shown, about 3- 10-4, 
is for a density of 1-10! states/cm?. For 2-10! states 1 cm?, one pre- 
dicts a ~ 6-10~*, which is reasonably close to the 4-10-* observed at 
this density.** Tentative results for other devices*! have yielded similar 
low-frequency a’s for interface-state densities up to 2-10" states/cm?. 
If true, a reexamination of the contribution of interface states would 
be in order. However, before this is attempted, it is essential to ascer- 
tain the interface-state density not at midgap but at V sg (see Appendix 
B) for the devices whose a’s are being measured. Lateral inhomo- 
geneities*® may also be contributing to incomplete transfer. 


4.4 The Conductively Connected Charge-Coupled Device (C4D) 


The C4D* is illustrated in Fig. 11. Diffused p-regions connect the 
depletion region under the capacitor plates as compared to the usual 
CCD arrangement of closely spaced capacitor plates. An ion-implanted 
n+ barrier is used to provide the directionality so that two-phase 
operation is achieved. From Fig. 11, 1t can be seen that transfer occurs 
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Fig. 11—Structure of the conductively connected, charge-coupled device (C4D). 


by a two-step process. The charge is stored under that portion of the 
capacitor plate which sees the n-silicon substrate (labeled C's in Fig. 
11), and the first step of the transfer, a CCD transfer process, is from 
this depletion region to the adjacent diffused p-region. The second step 
of the transfer is from the p-region over the n+ barrier to the next 
inversion region, a bucket-brigade transfer process. 

In practice the p-regions will have finite capacitance to the substrate 
as shown dotted in Fig. 11. This has two effects. First, the p-region 
capacitance is in parallel with the adjacent drain capacitance C's during 
transfer and charge must transfer back and forth from the p-region 
to Cs during clocking. Even if no charge is transferred over the barrier, 
a charge corresponding to the p-island capacitance multiplied by the 
voltage change across C's will flow to Cs, thus assuring that the inter- 
face states in the storage region become occupied during each cycle re- 
gardless of whether charge has been transferred or not. This has been 
referred to as an automatic ‘‘fat zero.’’?*4 However, the fat zero in- 
fluences only the first step of the two-step process, and, as previously 
discussed (see also Appendix A), it is primarily the second step which 
contributes to the incomplete transfer. 

The second effect of the p-region capacitance is that it adds to the 
intermediate capacitor Cg. Since for best operating performance we 
Wish to minimize C'z, it seems that best results will be achieved by 
making the p-island capacitance as small as possible. In fact, in the 
limit of negligible p-island capacitance with respect to barrier capaci- 
tance, the C4D will have the same performance as a two-phase ion- 
implanted CCD with zero electrode spacing (see Fig. 9), and this will 
represent its optimum performance capabilities. However, it is the 
p-island capacitance when the p-region is at its most positive voltage 
with respect to the substrate that is of importance, and this corre- 
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sponds to the largest value of its capacitance. In most cases this will 
not be negligible compared to the drain capacitance C's. If we consider 
the other extreme case, when p-island capacitance is large compared 
to Cs, then the C4D will perform like a simple bucket brigade but with 
the modulation terms multiplied by the ratio of p-island capacitance 
to depletion region capacitance plus p-island capacitance. Also, inter- 
face-state effects are similar to those in the ion-implanted two-phase 
CCD. Hence, C4D operating characteristics will lic between those of 
the simple bucket brigade and those of the ion-implanted two-phase 
CCD. 


4.5 The Stepped-Oxide and Tetrode Bucket Brigade 
The tetrode bucket brigade, first proposed by Sangster,?? is shown in 


lig. 12a. It was proposed?’ in order to reduce the drain conductance or 
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Fig. 12—(a) Structure of the tetrode bucket-brigade device. (b) Structure of the 
stepped-oxide bucket-brigade device. 
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feedback contribution to incomplete transfer, the effect known to be 
the dominant performance-limiting effect for bucket brigade at low 
frequencies. Sangster’ pointed out that because of the finite capaci- 
tance between the intermediate p-diffusion and the substrate, the 
improvement in performance was not as great as was hoped, but the 
drain conductance term ap was reduced by the ratio of the inter- 
mediate p-region capacitance to the storage capacitance. 

Recently,® one of the authors (CNB) has suggested a similar scheme 
to reduce ap referred to as the stepped-oxide bucket brigade and illus- 
trated in Fig. 12b. The gate of each IGFET is separated into two 
regions, the region nearer the source having a higher threshold than 
that near the drain. In this way the charge transfer process becomes 
a two-step process like that of the tetrode bucket brigade, but the 
intermediate capacitance Cz is made up only of the effective gate 
capacitance associated with the low-threshold region near the drain. 
In this way, the first step of the charge transfer is by a ‘“‘bucket- 
brigade”’ transfer process, and the second step is by a ‘“CCD”’ process 
for the stepped-oxide bucket brigade but by another ‘‘bucket-brigade”’ 
process for the tetrode bucket brigade. 

Like the other charge transfer devices treated here, the upper fre- 
quency limit for both the tetrode bucket brigade and the stepped-oxide 
bucket brigade will be limited by an intrinsic transfer rate term de- 
pendent on the channel lengths. Given the 10-micron assumption for 
all critical lengths, these two bucket-brigade registers will have the 
same upper frequency limit as the other CTD’s. At lower frequencies, 
the incomplete transfer parameter should become very small depend- 
ing on the ratio of Cg to Cp, and experimentally this has been found 
to be the case with measured values of a below 107*.* 


V. DISCUSSION AND CONCLUSIONS 


In the preceding section we discussed the coefficient of incomplete 
transfer, a, for CT'D’s in which the transfer of charge could be charac- 
terized as a single-step process or as a two-step process, and in which 
each step is either a “CCD” process or a ‘‘bucket-brigade’”’ process. 
All devices were found to have approximately the same high-frequency 
limitation dominated by the intrinsic transfer rate contribution to a 
under the assumption that they were all of the same geometrical size. 
However, the middle- and low-frequency performance of the single- 
step devices can be improved upon by using the more complicated 
two-step devices. One finds in two-step devices that the source is 
effectively isolated from the sink so that the primary contribution to 
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a comes from the modulation contributions associated with the second 
transfer. If now the intermediate capacitance is made much less than 
the storage capacitance, the modulation terms are reduced by approxi- 
mately the ratio of the intermediate capacitance to the storage capaci- 
tance. If this is achieved by making the storage capacitance large and 
keeping the intermediate capacitance fixed, better low-frequency per- 
formance will be accompanied by a lower high-frequency limit. 

It was also found that the bucket-brigade transfer process leads to a 
frequency-independent contribution to a at low frequencies. This re- 
sults from the requirement that charge must be injected from the 
diffused region into the channel of the IGFET. Conservation of current 
leads to a relation between the voltages of both sides of this barrier: 
V, in the diffused region and Vy at the beginning of the channel. For 
long transfer periods (low clock frequency) Va varies at 1/t and V> 
varies as log V4: hence, V, varies as log é, or essentially independent 
of time and of clock frequency. 

The analytic expressions we have obtained for the contributions to 
the coefficient of incomplete charge transfer a of both single-step and 
two-step transfer devices are quite general and provide a good qualita- 
tive and a reasonable quantitative prediction of the dependence of a 
upon the various device parameters. If more accurate quantitative 
results are desired, particular attention must be paid to channel-length 
modulation. This, however, is a static problem and should prove to be 
simpler than the dynamic problem of calculating a directly from a solu- 
tion of Q(é). 
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APPENDIX A 


Derivation of the Incomplete Transfer Parameter a 
A.1 Single-Step Process 


In a recent paper the authors" have shown that the incomplete 
transfer parameter a for a charge transfer device (ignoring for the 
moment the contribution of interface states) is made up of three terms 


a=a;+taptac 


where a; is the intrinsic transfer-rate contribution, ap is the drain-con- 
ductance contribution, and ag is the storage-capacitance modulation 
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term. For a single-step process, the intrinsic term is given by 


T Gx 
i — exp |] — oat (16) 
i i Cs 


where g, is the transconductance of the transfer mechanism joining 
C's to Cp (see Fig. 6) and 7 is the transfer time. As a first approximation 
for square-wave clocks, gm ~ dI/dVg and dt ~ (Cs/I)dVs, so eq. 
(16) becomes 

I 

a; ~ 17 

L, (17) 

where J, is the transfer current at the beginning of transfer. 
The drain conductance term, ap, is given by 


ee 
JrU s (18) 





ap = 
ImC p 


where g, 1s the reverse transconductance of the transfer mechanism, 
and both g, and gm are evaluated at time 7. Since in eq. (4) it is only 
the constant 8 which can be dependent on drain voltage Vp, we can 


write 
ol I / dg 1 dL. 
eH) EE) ow 
dV p B\dVp L.dVp 








where we have taken 8 = ZuC';/L. [see eq. (4) ]. 

According to the assumptions made here and listed in Table I, 8 
can only vary through channel-length modulation. However, the exact 
form of this variation has not been theoretically established. The most 
common assumption is that the channel length Lg is shortened by an 
amount equal to the one-dimensional depletion-layer width Lp asso- 
clated with the voltage on the IGFET drain diffusion. Even for an 
IGFET with source shorted to the substrate, this assumption is some- 
what crude, and for the bucket brigade with the source p-islands re- 
verse biased with respect to the substrate, it is expected to be even 
less accurate. For the CCD, an additional problem arises in that the 
transition fringing-field region under high lateral electric field connect- 
ing one charge storage region to the adjacent one extends under both 
capacitor plates so that less ‘‘channel-length modulation”’ will occur 
than for an equivalent bucket brigade.’ ”° This effect is particularly 
important when an ion-implanted barrier is used because, owing to 
the higher doping density in the barrier, channel length associated 


170 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1973 


with the barrier will change much less than the length of the adjacent 
storage well. An accurate treatment of this problem of calculating 
dLc¢/dV p is beyond the scope of this paper. However, in Appendix C 
we briefly outline an approximation which can be used. 

Returning to eq. (19), and inserting J from Eq. (4), we have 











B 2kT\ 1 dL. 
Jr = 5 Va (V. a —] L. iV, (20) 
Similarly, from eq. (4) we can derive gm to obtain 
ol KT\ dV a 
fe op, (vs ? —) aVs (21) 


lor a “‘bucket-brigade”’ transfer process, OV4/0V sis [1+ (kT /qV a4) ]7 
from eq. (11). It is also easy to show using eq. (11) that the effective 
C's for a bucket brigade-transfer process is (see Fig. 4) 


Va 


Cg = ——_§_— Cg: 22 
S !) + LT G ( ) 
Vat ee 
Hence ap Is given by 
Vit 2kT 
2kT gq 
Cy (V. = —) +CgVa 

q Vit kT 

: 1 dL, fas 
—— 2Cp Le dVp 


Given a calculation of V4 as a function of transfer time, the drain 
conductance contribution can be determined from eq. (23) by using 
the appropriate values of C, and Cg. For a CCD transfer process, Cp 
is zero, and in the limit as V4 tends to zero (i.e., at long transfer times 
or low clock-frequency operation) ap tends to zero linearly with V4. 
For a bucket-brigade transfer process, C’, is finite and ap will tend to 
a constant value in the limit of small Va. 
The capacitance modulation term, ac, has been shown! to be 


1 —*) 


Cs dQ. (24) 


ae = Qs( 


where Qs is the transferable charge remaining on C's at time 7. Usually, 
however, C's is made up of several contributions only one of which is 
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modulated by Q.. In this case if the particular capacitance is C's; with 
a transferable charge Qs; remaining on it, 


if 1 dCs; 
— Qs (— ~). (25) 
Cp \Cgi dVp 








In all the CTD’s of interest here it is the gate capacitance Cg which 
is modulated due to the channel-length modulation effect. Since this 
capacitance is linear in channel length, we obtain 


_ “= ( 1 =<) 


_ en 26 
ee Nas ) 








The one extra contribution to incomplete transfer which will be 
considered here is that due to interface states. It has been found! 
that such states lead to two terms, one due to interface-state capaci- 
tance modulation and the other due to the modulation of the dynamics 
of trapping and detrapping during transfer. For square-wave clocks 
and a large circulating charge or fat zero, it has been shown that the 
latter contribution will be small.*! Thus we will consider here only the 
capacitance modulation term ag.gg. If Qgs is the net charge trapped 
in interface states after transfer, then 





1 dCgg 
; (27) 


ac,ss = Qss(—— dQ 
SS o 


where C'gs is the interface-state capacitance. If we assume that this 
capacitance is also modulated by the channel-length modulation, then: 


Oss ( 1 ne 
es . 28 
ac,SS Cp io dV ( ) 





We have used this expression in a previous article.“ On the other hand, 
Tompsett”? has reported that variations in edge effects and in capture 
during transfer are important. These are included in (27); however, 
whereas the channel-width modulation leading to (28) can be greatly 
reduced by using two-step rather than single-step transfer processes 
(see below), edge effects, which depend on the initial amount of charge 
Q@., present, are nearly unaffected by merely increasing the number of 
transfer steps. Reduction of edge effects can come only by making 
the effective size of the well less dependent on Q,, e.g., by ion-im- 
planted barriers’ or diffused regions, or by keeping the carriers away 
from the surface, e.g., by storing charge in diffused islands or buried 
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channels. Changes in capture during transfer can be greatly reduced 
by using zero-gap technology.®? 

It should be noted that apart from the intrinsic contribution a;, the 
three contributions given above, ap, ac, and ac,ss, are all proportional 
to (1/Le)- (dLc/dV p). Hence, even though our estimates of dL¢/dV p 
tend to be crude (see Appendix C), the relative size of these contribu- 
tions can be ascertained more precisely. 


A.2 Two-Step Process 


In a two-step transfer process as shown in Fig. 6b, the charge first 
transfers from C's to an intermediate capacitor referred to as Cz, 
then to the drain capacitor Cp during a single transfer time 7. In this 
case transferable charge is left behind on both Cg and C's, and the 
derivation of a becomes somewhat more complex. The main advantage 
of using a two-step process is that all of the contributions to a due to 
channel-length modulation can be reduced since modulation of C's 
by the voltage Vp on Cp is a second-order effect (i.e., variations in V p 
modulate slightly the value of Vz, which in turn modulates to a much 
lesser degree the value of Vg). For this reason, we will ignore here the 
channel-length modulation effects for the first step of the two-step 
process. 

Referring to Fig. 6b for a definition of terms, the transferable charge 
on the first storage capacitor at time ¢ is 


Vs VSS1 
Qs = | CsdV + CsgidV. (29) 
V 


S,0 Vssi,0 


In this expression we shall be able to ignore the effects of channel- 
length modulation, these being of second order in a two-step process. 
However, channel-length modulation effects are important for the 
charge Qz on the intermediate capacitance C'z. We write 


VB Vss2 
Qe = | CpdV + CggodV (30) 
V Bo V ss2,0 


where Cgs2 is only the interface-state term subject to first-order 
channel-length modulation. We can solve for the effective a for the 
first step of the process a1 (see Appendix B) 


_ dQs 
- 


Qi; = exp |- | 2m at 
0 Cs 





a1 


= a1i + aic,ss (31) 


where 


TABLE [J]—Masor CONTRIBUTIONS TO INCOMPLETE CHARGE TRANSFER IN TERMS OF SMALL-SIGNAL 
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Modulation, ac,ss 
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TABLE [V—Masor ContTRIBUTIONS TO INCOMPLETE CHARGE TRANSFER IN TERMS OF APPROXIMATE 


Type 


Intrinsic, a; 


Drain Conductance, 
ap (Feedback) 


Capacitance 
Modulation, ec 


Interface State 
Capacitance 
Modulation, ac,ss 
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and aic,ss 18 given in (65). For the second step, however, the net cur- 
rent into Cz 1s the difference between J; and J, so that the equation 
to be solved is, using eq. (81), 


d ‘i © Om dV dV 
— = exp | - | 2 at | = gue + ga (32) 
dt C's 0 C's dQ)» dQ 


where 


— dQp 
ares 


(Details for arriving at an expression similar to (32) are given in 
Appendix B.) From eq. (30), 


(33) 





a2 

















dV 1 YB AC 1 Vssz dC 
ten a gee Save (34) 
dQ) o Cp Cx V Bo dQ « Cp Vss2,0 dQ. 
and | 
dVp 1 
= —- (35) 
dQ. Cp 
Solving for a2, in a similar manner to that used in Ref. 10, 
2 = dei + aep + aec + a2¢c,88 + oi (36) 
where 
"Qm2 TI 
Qe: = exp | — am at = — (37) 
| : Cx Ts, 
is the intrinsic transfer rate term, 
Jr2 Cp 
Aap ~ oo (38) 
° JY m2 Cp 


is the drain conductance or feedback term, which can be evaluated 
using eq. (23), 


1 dCpz Qe 1 <<) 
a cs Romicriry (eat 39 
Sea E wa Ci (- dV p Sa 


is the capacitance modulation term, 




















1 dC gg Qsso/ 1 dLee 
ase.ss ~ Qsae| l=) 
Css2 dQo Cp \Lez dVp 
is the interface-state capacitance modulation term, and 
/ mi C 
ai = a, (41) 


Jm2 C's 
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reflects the fact that charge transfers from C's to Cg in addition to the 
transfer of charge from Cz to Cp. Realizing that the total incomplete 
transfer parameter is a1 + a, the final result for a two-step transfer 
process can be written in terms of the parameters of the individual 
transfer steps as 


a = a2 + aap + arc + o2e¢c,ss + aic,ss + a1i (1 = in =). (42) 
J m2 C's 
Referring to eqs. (38), through (41), it is apparent that a significant 
improvement in a can be obtained in a two-step process by making 
Cp much smaller than Cy, provided the edge-effect contribution?’ to 
aic,ss can be reduced. 
The results of this appendix are summarized in Tables III and IV. 


APPENDIX B 
General Derivation of the Incomplete Transfer Parameter a 


In a previous paper" we outlined the details of the derivation of our 
general expression for a for the single-step process. In this appendix 
a similar derivation of a is given for the two-step process. The general 
derivation for an m-step process should then be straightforward to 
devise if needed. 

Referring to Fig. 6b, as the size of the initial, transferable charge 
@., 1s varied, both Qs(t) and Qz(t) will also vary. If Qn(t) = Q. 
— Qs(t) — Qz(t) is the transferable charge on the drain, then a for 
transfer process may be defined as 











d 
fea (43) 
dQ. 
It follows that 
a =a, tay (44) 
where 
OT = Ws (45) 
dQ. 
referring to J,, the current of the first transfer step, and where 
d 
ao > ae (46) 
dQ. 


which refers to both J; and I2, Z2 being the current of the second 
transfer step. As in the case of single-step transfer, we shall derive a 
differential equation which can then be solved for a. 


INCOMPLETE CHARGE TRANSFER V7 


As before," we first assume that we can write the currents governing 
the transfer of charge in the following form 


I, = I1(Vs, Vs, Ve, Vissi) (47) 
and 
I, = In(Vze, Vo, Ve, Vso) (48) 


where the additional voltages Vgg; (¢ = 1,2) are just the voltages 
induced by the trapped charges Qss;: on the effective capacitances 
Css: of the traps (Vc is the clock voltage). Qs, Qz, and Qp are related 
to Cs, Cz, Css: (¢ = 1,2), and Cp as follows 


Vs Vssi 
= | CsdV + CygsidV (49) 
Vis.o Vissi, 0 
VB Vss2 
Op = | CpdV + Csg2dV (50) 
VB,o V ss2, 0 
VD 
Op = | CpdV = Q. = Os = Op. (51) 
VpD,0 


From (49) to (51) and using (48) to (46) one can obtain the following 
relationships: 


dV 1 Ys ac 1 Yssi dC 
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dQ. Cz Cp Vzio dQ o Cz Cz V ssz,0 aQ)o 
dV 1—a, —- 1 Yo dC i 
ee se ae “= | at ome (54) 
dQ. Cp Co Jvp0 Qo Cp 
where we have made use of the following definition: 
dV gsi . 
agsi = Csgsi dQ, ) (¢ = 1,2). (55) 
A differential equation for a = a, + a, can now be derived at 
once. Since dQs/dt = — J, and dQz/dt = I, — Iz (conservation of 
charge), it follows that 
da dls dV z adVp OL» dV sso 
Se ee (56) 
dt dQ o dQ o dQ o OV sse dQ 


where as before gmz = O12/0Vz and grx = — O12/0V p. Inserting (53) 
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to (55) into (56), replacing a. by a — a, and solving for a, one finds 
that 


t 
a(t) = | ax(t’) E(t, t’) dt’ 
0 
t VB dC zp 
+| (| iv) E(t,t’)dt’ 
0 VB,0 dQ. 
t 
+ | (Gr2C B/G mC pd) H(t,t’) dt’ 
0 
rs I ( (1 ol, Cz 1 
0 ice OV sseC ss2 Gme 


Vss2 dC sso 
oe | wv) E(t, t')dt’ (57) 
V 


SS2,0 0 














where EH (¢,t’), the suppression factor, is defined by 


E(t,t’) = exp (-| au /ra(t")) [rll (58) 
and 


L/ro(t!) = gmat’)/Ca(@’). (59) 


The second, third, and fourth integrals in (57) can be evaluated at 
once assuming the factor multiplying E(¢,t’) varies more slowly than 
E(i,t’). This yields 


t VB 
a(t) = | a(t’) E(t,t’) dt’ + | dC 3/dQ dV +- gr2C B/Gm2C 
0 V B,0 











dl. Cz 1 a dC gs s2 
OV sse Csg2 Gme Vss2,0 Qo 


These terms are discussed in the text, in Appendix A, and in a previous 
work.!! 

To evaluate the first term in (57) we must calculate a, (¢). To do this 
we derive and then solve a differential equation for ai(¢). Since dQs/ 
dt = — J, it follows that 


day dl, dV g dV ep Oly dV gsi 
et ee, = = 4 —— A ° 
dt dQ dQ). dQ. OAVsgi dQ, 
Before inserting (52) and (53) into (61), we should note the relative 
magnitudes of the contributions to a1. dVz/dQ. is bounded by a/Cz 


-- aS s2 (1 av. (60) 











(61) 
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[see (60) ]. Hence the contribution to a, of the term g,1dV 2/dQ, is of 
the order of a-a;, which is quite negligible. As explained in Appendix 
A, variations in C’g are second order in AQ, and hence can be ignored 
in (52). In Ref. 9, we also pointed out that dC ss/dQ. was the dominant 
interface-state contribution to a, the term in agg being much smaller. 
Here dC ggi/dQo and dC'gs52./dQ, are the corresponding dominant trap- 
ping terms by the same reasoning. Hence, the third term in (52) as well 
as the fourth term in (60) can be dropped. This reduces (61) to 





d ie wie £ V8St cde 
aie a acon ed Say: (62) 
dt C's Cg Vssi,0 dQ) o 
the solution of which is 
at; (t) = a1;(t) + aic,ss(t) (63) 
where 
t 
cas(l) = exp (— | (oms/C send” (64) 
0 
and 


t t 
aic,ss(t) = | dt’ exp (- | (oms/C sent” 
0 t’ 


a Vssi dC 
(= | — av) (65) 
Cs Vssi,0 dQ t’ 














Vssi dC 
re | wr av (66) 
Vgsi,0 aQ 
1 dCgs1 
= 67 
ue Cssi1 dQ. oe 


Inserting (63) for a;(t) into (60) we obtain an explicit expression for 
a(t) as desired. 

Our expression for a(t) can be simplified somewhat. If aic,ss (66) 
can be assumed to be slowly varying near turnoff, then we can use the 
fact that fH dt = 1 [see (58) and (59) | to write 


t VB 
a(t) = | arni(t)E(U)dt? + | dC e/dQ AV + gr2C B/G mC dD 
0 VB,0 


Vss1 VSsse 
=F dC gs1/dQodV + dC gs2/dQ.dV. (68) 
V ssi, 0 V ssz2,0 
In (68), we have dropped terms proportional to ags1 and agg: as dis- 
cussed above. 
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To simplify the first term of (68) we must proceed more carefully. 
Using (64), we may write this intrinsic term as 


au(t | dt’ exp (- | dt" (9m2/C a — dma/Cs) )+ (malt) /Calt). 


This expression can be handled in various ways. In Appendix A we 
have assumed that gmo/Cz > gmi/Cs in which case one can multipy 
and divide the integrand by (gm2/Cs — Qmi/C's), assume (Gm2/C'p)° 
(gm2/Ca — Ymi/Cs)— is slowly varying, and integrate as before to 
obtain the intrinsic term of a 


1i° (1 — OmiC B/Gm2C gs) os oi(1 + GmiC B/Gm2C's). 


Alternatively, if Qmo/Ce = aof(t) and gmi/Cs = aif(t), then the 
integral may be performed without approximation, and one finds 


ani = - - {1 — |-f ie as) seat |} 


Finally, in the extreme case that gm2/Cs = 9mi/Cg one has for this term 





1; loge (1 /a1;) 


which varies typically from 5a; to 9a, for a; from 0.007 to 0.0001. 
Being the intrinsic term, it dominates in a only for relatively high 
clock frequencies." 


APPENDIX C 


In this appendix we consider several approximations which can be 
used to determine dLc/dV py. One approach is to simply define an 
appropriate one-dimensional depletion layer width Lp assumed to 
vary as the square root of voltage Vp, and assume that the effective 
channel length, LZ. is Lg — Lp. This approximation, valid only when 
the oxide thickness (times the ratio of the dielectric constant of the 
oxide to that of the semiconductor) is much less than the space-charge 
width, we shall use when a channel in the substrate material empties 
into a heavily doped diffused region, as is the case for the simple, 
tetrode, and stepped-oxide bucket brigade. Under these conditions 
we find using eq. (4) 


1 aL, Lp 1 
L.dVp 2(Le — Lp) |Vo| 


If, on the other hand, one has charge transport over a barrier into an in- 
verted region at the surface of the substrate material, the length modu- 





(69) 
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lation will be reduced approximately by a factor of (1 + N2/Ns5)4, 
where Nz is the doping in the barrier region and Nz is the doping in 
the substrate. This may be applied to the use of the C4D and the ion- 
implanted-barrier two-phase CCD. Thus one finds 

1 dL, Lp J 1 


EedVi. 2 Ge = hs). Veli Ne Ne 








(70) 


where Lp refers to the substrate material. Estimating dL./dVp by 
equating it to the modulation in the size of the depletion region should 
represent an upper limiting case. That is, the actual dL./dV p should 
be less than the value predicted by (69) or (70), making our resulting 
prediction for a conservatively larger than the actual result. 

For CCD structures in which charge transfers are directly between 
two inversion regions (e.g., as in the simple polyphase CCD and the 
stepped-oxide, two-phase CCD), the depletion-width approximation 
discussed above 1s clearly very unsatisfactory. We attempted to find a 
better estimate for dL./dV p in such cases but were unsuccessful. 

We hope that better approximations will be motivated by extensions 
of existing calculations. ”° 


Note: We point out that dl./dVp might be interpreted as the re- 
ciprocal of an electric field. The first field that comes to mind is the 
average (or peak) field at the silicon-silica interface between inversion 
regions. For the device dimensions used in examples in the text 
dL../dVp ~ (10° V/cm?)—!, whereas computer calculations?! suggest 
that the surface electric field is on the order of 104 to 10° V/cm. This 
disparity seems to rule out such a simple interpretation. 
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Quantizing Noise of AM/PCM Encoders 


By DAVID J. GDODMAN and LARRY J. GREENSTEIN 
(Manuscript received September 11, 1972) 


We consider the pulse-code-modulation encoder that contains a delta 
modulator for analog-to-digital conversion, and a finite impulse response 
digital filter that suppresses high-frequency components of the delta modu- 
latton signal. A PCM word generator produces fixed-length binary code 
words by rounding and amplitude limiting the filter output samples. The 
quantizing notse of the resulting PCM signal has four components: delta 
modulation slope overload noise, filtered delta modulation granular noise, 
amplitude overload noise, and word generator roundoff noise. We analyze 
the total quantizing notse for the case where the encoder input is a Gaussian 
random process and the digital filter impulse response is uniform (all 
coefficients equal). Such filters possess tmportant implementation ad- 
vantages and appear to be near optimal with respect to signal-to-noise 
performance. Our analysis resulis in curves which show the relationship 
of signal-to-noise ratio to filter order, delta modulation sampling rate, 
and PCM word length. 


I. INTRODUCTION 


A new approach to digital encoding of continuous waveforms 
employs digital hardware to unite the economy of single-integration 
delta modulation (AM) with the efficiency of pulse code modulation 
(PCM), A finite impulse response digital filter suppresses the granular 
noise component of the AM representation of a continuous signal, and 
a word generator truncates the binary coded filter output to produce 
PCM code words of desired length. This encoding methad controls the 
precision of the digital code hy means of the AM speed and the filter 
order rather than with the resolution of the multibit quantizer that 
appears In conventional PCM encoders, This is a desirable substitution 
in view of current technology in which the cost of high-speed digital 
circuitry is rapidly declining, 

This method of AM/PCM encoding, which was originally proposed 
by Goodman,’ has been applied to speech encoding by Freeny, et al.,? 

183 
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and to video by Kaneko and Ishiguro.* Previous theoretical results! 
focus on the filtering of the AM granular noise, but provide little 
insight into the important influence of AM slope overload and PCM 
amplitude overload on encoder design. Assuming the encoder input is 
a sample function of a Gaussian random process, the present paper 
analyzes the effects of the overload components of the quantizing 
distortion. It demonstrates that amplitude overload noise can be 
reduced if least significant bits of the filter output are truncated. 

We focus our attention on ‘uniform filter encoders,’”’ in which all 
filter impulse response coefficients are unity. Such encoders offer 
significant practical advantages, and they appear to be near optimal 
with respect to signal-to-noise performance. For such encoders, we 
show how performance varies with filter order, AM speed, and PCM 
word length, and we demonstrate the application of our results to the 
design of practical encoders. 


II, SIGNAL PROCESSING OPERATIONS 


The block diagram of Fig. 1 shows the operations involved in trans- 
forming the continuous signal y(¢) to a uniformly quantized M-bit 
PCM sequence. Digital logic may be added to convert this sequence 
to a nonuniform PCM format.® The single-integration delta modulator 
of Fig. 2 converts y(¢) to a sequence of pulses with amplitude +1 or —1 
at the rate f; = 1/r per second. The feedback loop is an ideal inte- 
grator with gain factor 6, while the up/down counter obtains a digital 
replica of x(t), the AM approximation signal. The output of the Nth- 
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Fig. 1—Encoder block diagram. 
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Fig. 2—Delta modulator. 


order digital filter is a weighted sum of consecutive AM approximations 
to y(t). 

Although the filter outputs range over a discrete set, the number of 
possible filter outputs is unlimited because there is no fundamental 
restriction on the amplitude range represented by a delta modulator. 
It follows that, with the PCM word length prespecified, an additional 
quantizing operation is required. This quantization 1s performed by 
the word generator which restricts to 2” the number of possible coder 
output words. The word generator introduces amplitude overload and 
it may also add to the granular quantization error by rounding off 
least significant bits of the filter output. 

The filter and word generator are controlled by a clock which causes 
coder output words to be generated at the rate 2W per second, where 
W is the bandwidth of the analog input. Hence, the data rate of the 
coder is 2/W bits/second and the PCM sequence may be decoded as 
if it were produced by a conventional encoder consisting of a 1/2W- 
second sampler and a uniform quantizer with 2” output levels. 


III. THE PCM QUANTIZATION LEVELS 


With the filter coefficients, g;, integers as in a practical realization, 
the filter output at ¢ = kr is the integer 
N-1 
Un = Ll Giwei (1) 
im0 
where {w;} is the sequence of counter outputs. Because w; = wj-1 + 1, 
the parity of the filter input alternates between even and odd at each 
AM sampling instant. It follows that if f,/2W, the ratio of AM sam- 
pling rate to PCM sampling rate, is an even integer, the parities of 
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Wk, Wk-1ly °° *, Wk-n41 are invariant at the PCM sampling instants. 
Hence, the parity of wu; is the same at all PCM sampling instants. 
Because odd-parity filter outputs lead to an easily implemented word 
generator, we restrict our attention to encoders in which w; and ux 
are both odd integers at the PCM sampling instants. The filter coeffi- 
cients of these encoders satisfy conditions, derived in the Appendix, 
which do not severely restrict the set of available filter transfer func- 
tions. The conditions do, however, preclude uniform filter encoders of 
orders 4, 8, 12, ete. | 

If to is the encoder delay, the odd integer u; is a scaled approximation 
to y(kr — to). To determine the scaling factor, we observe that x(kr), 
the AM approximation to y(kr), is related to wz, by xw(kr) = dws. 
Further, since the filter provides relatively distortionless gain over the 
signal bandwidth, it expands the amplitude scale of w,; by approxi- 
mately the amount of the de gain, 

N~1 
T= 2b) 4. (2) 
i=0 

Thus, (6/J)u;, is an approximation to y(k7 — to) and, with u, ranging 
over odd integers, the signal levels represented by the input to the 
word generator are in the set 


6 6 
7 Bart (3) | 
with quantizing step size 26/J. Because the scaling by 6/J is approxi- 
mate, we admit an additional scale factor, y, which brings the PCM 
representation optimally close to y(t) in the mean square sense. The 
actual step size of the filter output is therefore 
26 

dy = rae (4) 
In Section 7.5, we show that y, which depends on g;, is close to unity 
for encoders of practical interest. 

Figure 3 shows the mapping of the filter output into M-bit code 
words. To eliminate a information bits from the binary representa- 
tion of u,, the word generator truncates the a + 1 least significant 
bits. (With u, odd at PCM sampling instants, the least significant bit 
always has value one and hence conveys no information.) In the 
absence of amplitude overload, 


| wa, | < QMt+a _ | 
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Fig. 3—Word generator roundoff procedure. 


and the PCM code word consists of the M — 1 least significant bits 
of the truncated binary representation’ of u,; and the sign bit. When 


lun| > 2@t+e — 1, 


the transmitted code word is either the most positive or most negative 
M-bit word. A decoder recovers the integer code words of Fig. 3 by 
appending a one and a zeros to the least significant end of the PCM 
word. | 

With the truncation of each information bit, the step size increases 
by a factor 2 so that, with a information bits truncated, the PCM code 


T In Section V, we point out an advantage of the twos complement binary format. 
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words represent signal levels in the set 
d. d. do do 
mo ia = oem aoe ieee aida 1D (5) 
2 2 2 2 


with step size 


26 | 
d = 2d, = 2%*—v (6) 
I 
and maximum amplitude (2% — 1)2%6y/I. 


IV. UNIFORM DIGITAL FILTERS 


The value of J, the de gain of the digital filter, is crucial in determin- 
ing the character of the overall PCM quantizing noise. With the filter 
coefficients all integers, J may be regarded as a measure of coefficient 
quantization. A large value of J corresponds to fine quantization be- 
cause it allows considerable freedom in choosing g;. To obtain a filter 
transfer function that approximates with arbitrary accuracy the opti- 
mum transfer function with respect to granular noise,! an arbitrarily 
high value of J is required. On the other hand, amplitude overload 
noise increases rapidly with J because the dynamic range of the 
encoder is nearly proportional to [-1. 

The rapid increase in amplitude overload noise as a function of I 
leads us to focus our attention on the uniform filter, 


go = 1; 4 =O Ly eh IL, (7) 


for which J = N, resulting in the greatest dynamic range attainable 
with an Nth order filter with all coefficients of the same polarity. (We 
exclude from consideration filters with g; = 0 for one or more 7.) 
Reference 1 suggests that, for high sampling rates, the coefficients of 
the optimum filter with respect to granular noise are nearly equal and 
that the difference in granular noise rejection between this optimum 
filter and the uniform filter is marginal. This observation suggests that 
encoders with uniform filters, because they minimize amplitude over- 
load noise and produce near minimal granular noise, are nearly optimal 
with respect to total quantizing noise. Further support for this specu- 
lation is given later. 

In the frequency domain, the uniform filter transfer function is 
sin (rN f/f;)/sin (rf/f;) and the filter rejects increasing amounts of 
AM granular noise as N increases. So long as f;,/N 1s large relative to 
2W, the signal component of x(t) is undistorted by the filter; but, as 
f./N approaches 2W, distortion of in-band signal components becomes 
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significant, and overall performance deteriorates with increasing NV. 
Thus, if the advantages of very high-order filtering are sought, designs 
more sophisticated than eq. (7) are required. 


V. IMPLEMENTATION 


Besides possessing noise-rejection properties, uniform filters admit 
considerable hardware economies relative to other designs. With all 
coefficients unity, no multiplication is required, and each filter output 
is merely the sum of N successive counter levels. Therefore, one may 
implement the uniform filter as a resettable accumulator, thereby 
eliminating the delay line of Fig. 1, as well as the multipliers. To obtain 
a PCM sample, the coder sets the accumulator to the current level of 
the up/down counter and adds to the accumulator the next N — 1 
counter levels. 

Because the addition of N numbers is required only once for each 
PCM sample, and because f,/2W, the number of AM samples per PCM 
sample, is generally much greater than N, it is possible to time-share 
a single accumulator among many signal channels. With inputs pre- 
sented to the accumulator at the AM rate, the number of channels 
sharing a single accumulator may be as high as f,/2WN. Hence, in 
terms of hardware requirements, the filter order, NV, determines time- 
sharing capacity rather than the number of circuits necessary to 
realize a single encoder. 

In addition to adding counter levels into an accumulator and 
truncating least significant bits of the sum, the encoder must detect 
amplitude overload and generate the most positive or most negative 
code word when the word generator is overloaded. 

It must also restore the proportionality of the counter level, wz, 
to the AM approximation, x(kr), after each instance of counter over- 
load. The wrap-around property of twos complement arithmetic 
ensures this proportionality whenever |x(kr)| < (2%+« — 1)6. On the 
other hand, a saturating counter would require special measures to 
restore tracking after each instance of counter overload. 


VI. ENCODER PERFORMANCE 


6.1 Figure of Merit and Design Specifications 


An ideal decoder of the encoder output sequence obtains @(t) ( defined 
in Section 7.1), a delayed, noisy approximation to the analog input 
y(t). We define the quantizing noise power of 7(t) to be 


Nr = Etlyg® — yt — to) P} (8) 


190 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1973 


. i 
GRANULAR | OVERLOAD 
J 


SLOPE = 1 
(10 dB / DECADE) 


SIGNAL—TO-—NOISE RATIO IN DECIBELS 





(7/8)? (18)? ont 
(o/8)* LOG SCALE 


Fig, 4—Relationship of design objectives (7n,7 Rr) to figure of merit (2). 


where Ff denotes expectation, to is the encoder delay, and y is the 
mean-square optimum scaling factor. For each digital filter, PCM word 
length, and AM sampling rate there is a unique combination of values 
of 6, the AM step size, and a, the word generator parameter, that 
results in minimal Ny. We choose as a figure of encoder merit the ratio 
of signal power to this minimum noise power,! 


| min Ny 


a,d 


(9) 


In the design of a practical encoder, typical specifications include a 
signal-to-noise ratio design goal, Zp, and a range, Zr, of input powers 
over which the actual signal-to-noise ratio must equal or exceed Zp. 
The practical significance of our figure of merit is found in the 
approximation 


4x Lpn+ Ze (10) 


' Nr is a convex function of 6 and a In our numerical work we have used simple 
search techniques to find ming,s Nr. 
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where each quantity is measured in decibels. For example, an encoder 
for which Z = 55 dB will actually attain this signal-to-noise ratio for 
a single level of input power, and will maintain a signal-to-noise ratio 
of 35 dB or better over a range of 20 dB in signal power. 

Equation (10) is derived from Fig. 4, an approximation to the 
dependence of signal-to-noise ratio on input level. If (¢/6) exceeds 
(c/5)opt, the optimum ratio of rms input to AM step size, overload 
noise predominates in the distortion and the signal-to-noise ratio falls 
rapidly as o? increases. On the other hand, with (0/5) < (o/8)opt, 
granular noise predominates and, with 6 fixed, is essentially indepen- 
dent of c. Hence, the signal-to-noise ratio 1s proportional to o? in the 
granular region. 

If the ensemble of input power levels is log-normally distributed, as 
in models used for speech signals,* 10 log (¢/6)? is a normal random 
variable, the mean value of which we denote by 10 log (¢/6). Hence, 
the probability that the signal-to-noise ratio exceeds Zp 1s maximum 
when dp, the design value of the step size, 1s chosen such that 


10 log (¢/5p)? = 10 log (¢/6)% pt — $Zp. (11) 


That is, 10 log (¢/6p)? is the midpoint of the design range of length Ze. 


6.2 Performance Characteristics 


Figure 5 shows a typical set of performance curves, computed 
according to the theory presented in Section VII. The curves pertain 
to 11-bit encoding of Gaussian signals having a truncated first-order 
Butterworth power spectrum, where the ratio of 3-dB frequency to 
cutoff frequency is 0.25. This type of process has been used to model 
band-limited speech.’ The performance curves show the figure of 
merit, Z, of uniform filter encoders of various orders as a function of 
f:/2W, the AM sampling rate expressed as a multiple of the PCM rate. 

The choice of a specific encoder configuration represents a com- 
promise between the advantages of low AM speed and low filter order. 
The nature of this compromise is illustrated in Fig. 6, which shows 
combinations of AM speed and filter order that satisfy two quality 
objectives. The broken curves relate AM speed to the maximum filter 
order consistent with sharing the accumulator described in Section V 
among 24, 48, and 96 signal channels, respectively. All design points to 
the right of a broken line are permissible for the given number of multi- 
plexed channels. 
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Fig. 5—Performance characteristics, M = 11 bits, Butterworth spectrum. 


6.3 Dependence of Performance on Design Parameters 


Figure 5 demonstrates two types of variation of Z with f,: Z rising 
with a slope of 20 dB/decade, and Z flat or decreasing slowly with f,. 
The first type of behavior occurs when amplitude overload is negligible 
and slope overload controls the optimum AM step size. In this case, 
the optimum step size varies approximately as 1/f,, and the decrease 
continues until amplitude overload becomes significant. When ampli- 
tude overload is the predominant overload noise, the optimum step 
size is essentially constant and the slightly negative slope of Z indicates 
that an increase in f, results in an increase in the granular noise corre- 
lation from sample to sample, leading to a greater proportion of the 
AM granular noise power in the passband of the filter. 

The flat portions of the curves represent transition regions to higher 


AmM/PCM ENCODERS 193 


FILTER ORDER N 
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Fig. 6—Design alternatives derived from Fig. 5, M = 11 bits, Butterworth 
spectrum. 


values of a, at which the increased roundoff noise of the word generator 
is offset by improved immunity to amplitude overload. As f, increases 
indefinitely, Z approaches the maximum signal-to-noise ratio associated 
with uniform PCM encoding of Gaussian signals. 

In Figs. 7 and 8, we see that the shapes of the characteristic curves 
are essentially invariant with the number of bits in the PCM code. 
In Fig. 7, which pertains to 11-bit encoding of signals with a flat 
spectrum, amplitude overload effects occur at points that are approxi- 
mately 10 log (2/28) = 18 dB higher in signal-to-noise ratio and 
further to the right by the factor 2!!/2° in sampling rate, relative to the 
corresponding points in Fig. 8, which pertains to 8-bit encoding of the 
same input process. 

Figures 5, 7, and 8 also demonstrate the effect of filter order. When 
f. is quite low and amplitude overload is negligible, Z increases mono- 
tonically with N. However, the value of f,; at which amplitude overload 
becomes significant decreases as N increases, and the earlier transitions 
from a=0 toa=1,a=1 to a = 2, etc., lead to the crossovers. 

Figures 5 and 7 relate to the same PCM word length but different 
signal spectra. The principal difference between the two sets of curves 
is a scale change of the horizontal axis. In Fig. 5, the axis is shifted to 
the left relative to Fig. 7 by the factor 1.6, which is the ratio by which 
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Fig. 7—Performance characteristics, M = 11 bits, uniform spectrum. 


the rms slope of signals having the uniform spectrum exceeds the rms 
slope of signals having the Butterworth spectrum. 


VII. QUANTIZING NOISE ANALYSIS 
7.1 Noise Components 


To reconstruct an analog signal from the sequence of word generator 
outputs, we first recover one of the integers shown in the ath column 
of Fig. 3 by appending a one and a zeros to the least significant end of 
each code word. We next multiply the sequence of integers by the 
nominal scale factor 6/J and denote the resulting sequence by {9;}. 
Finally, we perform ideal interpolation of { 9;} to obtain the continuous 


waveform 
: - sin 2rW(t — 7/2W) 
1) = 2. 4; ——— 


: (12) 
j=—oo 27 W (t P= 3/2W) 
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Fig. 8—Performance characteristics, M = 8 bits, uniform spectrum. 


Our purpose, in this section, is to investigate the difference between 
g(t) and the encoder input, y(¢), when this input is a sample function 
of a zero-mean stationary Gaussian random process. Our measure of 
distortion is the total quantizing noise power, N77, defined in eq. (8). 
Because y(t) 1s stationary, Nr is independent of time, and in the sequel 
we omit time arguments and subscripts from the notation of signals 
when there is no risk of ambiguity. Although y in eq. (8) is a compli- 
cated function of signal statistics and encoder design parameters, the 
introduction of the preliminary scaling factor, 6/J, leads to y = 1 in 
situations of greatest interest. (See Section 7.5.) The other constant 
in eq. (8) is to, and we observe that if the filter coefficients have even 
symmetry (g: = gn-1-:), the filter delay is 


to = (N —1)7/2, (13) 


one-half the filter memory span. 
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To derive Nv as a function of the encoder design parameters, we 
recognize 7 as the sum of a signal term and four noise terms. We begin 
by writing 


6 
y a T (u 2 Nw) (14) 


where wu is the filter output and n» represents the difference between 
the input and output of the word generator. In studying (6/J)u, the 
filtered AM signal, it is customary to identify granular and slope over- 
load components of the AM quantizing noise, in the manner indicated 
in Fig. 2. If we rewrite eq. (1) as u = g * (#/6), with * denoting con- 
volution and x denoting the approximation signal in the delta modu- 
lator feedback loop, we obtain 


6 
=u = oer =a (y + no + ns) (15) 
where ng and ns are AM granular and slope overload noise, respectively. 
The distortion introduced by the word generator, (6/J)nwy, may 
itself be resolved into two components, namely, n4, which accounts for 
amplitude overload, and nr, which represents the roundoff effect. 
These observations lead us to a representation of the quantizing noise 
signal as the sum of four noise components: 


(yg — y) = [B+ + no) — y+ Bens + ama + ame. (16) 


The term in square brackets is the filtered granular noise, and the 
remaining terms are slope overload noise, amplitude overload noise, and 
word generator roundoff noise, respectively. In this paper we evaluate 
the expected square of eq. (16) by assuming that the expected product 
of each pair of terms is negligible relative to the sum of the two mean 
square values. Thus we express the total quantizing noise as the sum 
of four noise powers 


Nr=NetYNgtyNat Ne (17) 


in which each term is the expected square of a term in eq. (16). 

When the total quantizing noise is low, we are justified in approxi- 
mating the average cross products of eq. (16) by zero because: (2) 
granular noise and roundoff noise are zero during overload bursts; (72) 
each type of overload occurs with low probability and the probability 
of their joint occurrence is negligible; and (zi7) we have found that 
|E(ngnr)| is many orders of magnitude lower than Ng + Nz when 
M 2 4. 
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7.2 AM Granular Noise and Slope Overload Noise 


Expanding the square of the term in brackets in eq. (16), we obtain 


y NaI 
Ne = Ey’) - os Dd, giltey(tr — to) 
i=0 
y\2 Voi No 
. G) LL gigies(tr — jr) (18) 
7=0 7=0 


in which the AM correlation functions, 
Riy(7) = Ely@O2@4+7)], Rar) = Ela@rc + 7)] 


are derived in Ref. 8, under the assumption that overload never occurs. 
To use the results of Ref. 8 and also account for overload, we should 
multiply Neg in eq. (18) by the probability that overload is absent. For 
the applications that interest us, the probability is greater than 0.99, 
and we adopt eq. (18) as an approximation to N¢ that overestimates 
the granular noise component of Nr by no more than one percent. In 
Section 7.4, we similarly overestimate Nr. 

To compute Ns, we adopt the assumption of previous authors’? that 
essentially all of the AM slope overload noise power is in the signal 
band of y(t) so that Ns = E{[(g/I) * ng ]?} = E(ng). In our numerical 
analysis, we have followed Protonotarios,? who derives E(ng) as a 


function of 
5 dy} —4 : 
s=-\2|=1) ; (19) 
T dt 


the ratio of the maximum slope of x(t) to the rms slope of y(¢). For 
high values of S, Ns is proportional to S~> exp [—38S? ]. 
7.3 Amplitude Overload Noise 


The maximum output of the word generator is 2¢(2” — 1). Assum- 
ing there is no granular or slope overload noise during amplitude 
overload intervals, 


) 
nam aay — 5292" — 1); g*ey > 2212” — 1)6 


6 
= Fey +2202" — 1); gt <= 2°" — 1)0 
205 lg «y| S222" — 1)8. (20) 


Because (g/J) * y is a sample function of a zero-mean Gaussian process 
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with variance 


; | N-1 N-1 
or =— dL Dd gig;)Ry,(ir — Jr), (21) 


ee t=0 j=0 


the mean square value of eq. (20) is 


Na 


2/n) | o p(x — A)? exp (—42?)dx 


A 


or{(1 + A) erfe (2-14) — (2/m)3A exp (—3.42)} (22) 


I 


in which A is the amplitude overload factor, 
2°(2™@ — 1)6 


lor 


A= (28) 
7.4 Word Generator Roundoff Noise 
With wu the filter output and z(:) the mapping shown in Fig. 3, 
62 
Nr = pee oe) al (24) 


Because u is an odd integer, we have the binary number representation 
of u > 0, 


u>= >, b;2* -+- 1 
t=] 
where 6; = 0 or 1. The word generator truncates baba_1 --- b11 from 


this representation and z is obtained by replacing these digits with 
10--- 0 = 2@. 


Hence 
a > 6,2? + 2° (25) 
t=at+l1 
and 
z2—u=Q; a = 0 
= 2*—] — Do b,2'; wel. (26) 
i=l 

For u <0 the odd symmetry of Fig. 3 implies z(u) = — 2(—u). 


Hence z — wu is an odd integer in [— (2% — 1), 2* — 1] and the ex- 
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pectation in eq. (24) is a weighted average of the integers 


1?, 37, --- , (2% — 1)?. From this observation, we immediately obtain 
the bounds 

6? 6? 

eee TOE: 2 a 2 

7 SNe rR (2 1)? (27) 


With a = 1, the bounds are equal and Nz = &/T’. 

For a > 1, we evaluate Ne only for coders with uniform digital 
filters. All odd integers are possible outputs of such filters. That is, 
Pr{u = 2n + 1} > 0 for all n and, for low-noise encoding, this prob- 
ability is quite nearly constant over a set of 2°~! consecutive integers. 
When scaled to the amplitude range of y(¢), such a set of filter outputs 
lies in an interval of length 


26 
7 2e-ly = d/2, (28) 


a small fraction of o (typically of the order of 4¢/2”). For M2 5, 
the envelope of Pr{u = 2n +1} has approximately the Gaussian 
shape of the probability density of y(t) and a piecewise constant ap- 
proximation to this density over intervals of length d or less leads to 
highly accurate expressions for quantizing noise power.!° 

Over intervals of length d/2, (2 — u)? takes on all the values 1?, 
37, +++ (2% — 1)? either in ascending or descending order. Hence, the 
piecewise constant approximation to Pr{u = 2n + 1} reduces the ex- 
pectation in eq. (24) to an unweighted average of these 2° integers, 





Np = ogre > (21 — 1) de Solio (29) 
= — JO Na 1— = — ° 
“7? = I 3 


Noting that eq. (6) admits the expression 


: (=) _ qd? 
TNT) > 4(49)’ 


we summarize the results of this section as follows: 


Nr = 0; = 0 
qd? 

~ 16’ a= l 
qd? 


i (1 — 4-2); a = 0, uniform filter encoders. (30) 
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| The last line indicates that Nr—- @/12 asa—«. This limit is the 
granular noise associated with instantaneous PCM encoding of samples 
ranging over the continuum. 


7.5 The Optimum Scale Factor 


To complete our quantizing noise analysis and establish the validity 
of prescaling the word generator output by 6/Z, we show that y, the 
additional scaling factor that brings the amplitude of §(¢) optimally 
close to that of y(t) in the mean square sense, is nearly unity in designs 
of practical interest. Specifically, we derive the inequalities, 


(1—b)yo<¥<Y0 (31) 
where b is of the order of magnitude of the noise-to-signal ratio and 


NV=1l N-1 


2 a £igihyy(tr —_ to) 


t=0 j=0 
0= 
Y N 


(32) 


—1 N=1 
De £igjlyy(ir — jr) 
t=0 j=0 

| Clearly as f, >, all of the covariances in eq. (82) approach o? and 
hence yo-> 1. In all of the numerical examples considered in Section 
VI, the sampling rates are so high that yo is quite nearly unity. 
For example, all of the points plotted in Fig. 6 correspond to 
0.99 < yo < 1.01. 

By definition, y is the value of c that minimizes 


; 8 ° 
EL(cg —y)? J =F Ie «(y+ ne) —y + ens] | (33) 
The expression in square brackets on the right side is identical to eq. 
(16) with c replacing y and cno replacing the last three terms. Because 


we have assumed the correlation of mo and [ (cg/I) * (y + ne) — y] 
to be zero, we may rewrite 


El(c9 — y)*] = E E ae v| 4 cB(n2). (34) 


| Differentiating this equation with respect to c, and equating to y the 
value of c that causes the derivative to be zero, we obtain 
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Bly 2 + y+ no) | 
Be *(y + ne) | 4+ E(ni) 


] V-1 


F DL gilvsy(tr — to) 
i=0 
= N-1 N-1 (35) 
72 > a gig)Rec(tr = qT) +NgtNa+Nr 
t=0 j=0 


Reference 8 demonstrates that, for low-noise encoding, the approxi- 
mations R., ~ Ry, and R,z ~ Ry, + Ree [where R.-e(-) is the auto- 
correlation function of the unfiltered granular quantizing noise | are ex- 
tremely precise. Thus eq. (85) becomes 


1 N-1 
7 D. BR gir = to) 
40 


a ET 7 aaa 
7 DL 2 £igiltyy(ir — jr) + Nae + Ng t+ Nat+Ne 
im0 j=0 
where 
; N-1 N-1 
Ne=—D > 2:2; Reeltr Ih) 


[? i=0 j=0 


and is of the order of magnitude of Ng.1 We now recognize that the 
first term in the denominator of eq. (36) is ¢7[eq. (21) ], the power of 
the filtered signal component of the AM approximation. If we divide 
numerator and denominator by this term (and substitute >~; g; for J) 
we obtain 





Yo 
= 37 
Cor 5 (37) 
where we have defined 
pa A Nae Ne Na. (38) 


which is the order of magnitude of N7/o?. Equation (81) follows 
immediately from eq. (37). 
VIII. SUMMARY AND CONCLUSION 


Section VII presents the analytical steps enabling us to compute 
_ the figure of merit, Z, of a AM/PCM encoder. The rationale for using 
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this figure of merit is presented in Section VI (Fig. 4) along with 
results for some special cases of interest (Figs. 5, 7, and 8). The ultimate 
utility of these results is that they enable the designer to determine 
tradeoffs between AM sampling speed and digital filter order for 
specified values of encoder quality (e.g., Fig. 6). 

We should reiterate the conditions assumed for the encoder in deriv- 
ing our results. First, we have assumed that the digital filter output is 
at odd parity at every PCM sampling instant. Aside from simplifying 
the roundoff noise analysis, this condition appears to correspond to the 
simplest possible implementation of the PCM word generator. The 
primary design constraint it imposes is the prohibition of digital filters 
of orders 4, 8, 12, etc. 

Second, we have assumed that the digital filter has uniform co- 
efficients. This condition makes a complete noise analysis relatively 
straightforward and also leads to a simple filter implementation. 
I'urthermore, it corresponds to a robust design that appears to be 
near optimal in all cases of practical interest. Attempts to demonstrate 
the latter point quantitatively have foundered on the difficulty of 
assessing the roundoff noise power (Nr) when the filter coefficients are 
nonuniform. If we assume that, for any step size, d, the roundoff noise 
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Fig. 9—Uniform filter assessment. 


AmM/PCM ENCODERS 203 


power is nondecreasing with the number of bits truncated, we can 
derive an upper bound on Z. A comparison of this upper bound, Z, 
with the figure of merit resulting from uniform coefficients is shown in 
Fig. 9 for one case. Comparable results can be expected for other 
combinations of signal spectrum, PCM word length, and filter order. 


APPENDIX 
Odd-Parity Filter Outputs 


With the encoder delay equal to one-half the filter memory span 
Leq. (13) ], we consider symmetric coefficient sets 


gi = 9nN-1-1, 7=0,1,---,3N-—1; N even 


(39) 
gi = 9n-1-i; 7=0,1,---,#(N — 3); N odd. 


Such coefficients give equal weight to counter levels that are equally 
advanced or retarded with respect to y(kr — to), the input value 
estimated at ¢ = kr. Equation (39), when combined with eq. (1), 
implies 


}N-1 
Ut = ~ g;[ Wei + Wi—(N—1) +41} N even (40) 
inl 
$(V-3) 
Uk = SN We-yiw-ay HD iL We-i + We-w—-1) 4:1]; Nodd. (41) 
t=0 


With the counter levels, w;, alternating in parity, the two counter 
levels in square brackets in eq. (40) are of opposite parity because the 
difference in subscripts, N — 1 — 22, is an odd number. Hence, their 
sum is odd. On the other hand, the two corresponding counter levels in 
eq. (41) have the same parity, and thus an even sum, because 
N — 1— 2 is even with N odd. These observations lead to the 
following necessary and sufficient conditions for uz, ranging over the 
set of odd integers: 


Condition A: With N even, uz is odd if and only if there is an odd 
number of odd coefficients in the set go, g1, °°: , G4n-1- 

Condition B: With N odd, u; is odd at a PCM sampling instant if 
and only if gi~v_-1) is odd and the low-speed clock is synchronized so 
that PCM sampling instants occur when w;z-4:w-1) is odd, This syn- 
chronization can be achieved if the ratio of AM sampling rate to PCM 
sampling rate is an even integer. 
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In uniform filter encoders, Condition A is always satisfied when 
N = 2, 6, 10, etc. It can never be satisfied with N = 4, 8, 12, etc. For 
Condition B to be satisfied, the encoder must be synchronized such 
that w,z_v+1 (the first term entering the accumulator described in 
Section V) is odd when N = 1, 5, 9, etc; wz_w3.1 must be even when 
N = 3,7, 11, ete. 
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Scattering Losses Caused by the Support 
Structure of an Uncladded Fiber 


By D. MARCUSE 
(Manuscript received September 5, 1972) 


We consider an uncladded dielectric waveguide core that 1s held by 
dielectric supports. The radiation losses caused by the support structure 
are being considered. The analysis is sumplified by using a slab waveguide 
model held by slab-shaped supports. Only order-of-magnitude estimates 
are attempted. The radiation losses can be reduced by reducing the refractive 
index contrast between the supports and the surrounding medium with 
the help of an index matching liquid. The radiation losses remain large 
unless the index match ts sufficiently close. 


I. INTRODUCTION 


A typical optical fiber consists of a core, the refractive index of which 
is larger than that of the cladding material surrounding the core.! The 
cladding serves the purpose of keeping any outside influence, such as 
dust, at a safe distance from the core. The requirement of a lower 
refractive index for the cladding makes it difficult to find suitable 
cladding materials for one of the most promising core materials—fused 
silica. Fused silica is particularly useful as fiber core material because 
of its inherently low absorption loss. However, the refractive index of 
fused silica n = 1.46, is lower than that of most other glasses. In 
particular, there are as yet no low-loss materials suitable for fiber 
claddings, the refractive indices of which are lower than that of fused 
silica. The few available materials have absorption losses that rule out 
their use as claddings for low-loss optical fiber waveguides. 

It appears natural to ask whether a fiber without cladding could not 
be made. In principle a dielectric fiber waveguide works Just as well 
without a cladding if it could be suspended in air or in vacuum. But 
since no method of levitating the fused silica fiber without mechanical 
supports has yet been devised, the necessity exists of holding the un- 
cladded core by some kind of supporting structure. If the supports of 
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the fiber core have a refractive index lower than that of the naked core, 
the amount of light scattered out of the core by contact with the sup- 
ports may be tolerable. Since the supports touch the naked core only 
occasionally, their dielectric losses need not be exceptionally small. 
Solid materials of lower index than fused silica, but with much higher 
losses, do exist. 

The scattering losses of the supports could be reduced by submerging 
the core and the support structure in a low-loss index matching liquid.’ 
If it were possible to match the refractive index of the supports per- 
fectly, the submersion technique could eliminate all scattering losses 
from the supports so that only the losses of the index matching liquid 
and the heat loss of the fiber supports would remain. The heat losses 
of the supports would equal their average value averaged over the 
entire length of the guide. Because of the low filling factor, this average 
loss could be tolerably small. However, even this index matching 
technique encounters certain problems. It is unlikely that a low-loss 
liquid could be found that matches the refractive index of the support 
structure perfectly. But even if this were possible, the index match 
would work only at one fixed temperature so that scattering losses 
would still occur if the ambient temperature drifts from the design 
value. 

In order to explore the requirements for low scattering losses of a 
partially index-matched support structure, a model calculation is 
carried out in this paper. For simplicity we limit the discussion to the 
TE modes of the slab waveguide in the hope of obtaining order-of- 
magnitude estimates. The model to be investigated is shown in Fig. 1. 
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Fig. 1—Dielectric slab waveguide core of refractive index n, supported by two 
slabs of index n3. 
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The dielectric slab of refractive index n; is embedded in a medium of 
refractive index no. This surrounding medium can be thought of either 
as air or as a suitable index matching liquid. The support structure is 
simulated by two dielectric slabs attached at right angles to the 
waveguide core. The refractive index of these support pieces is 73. The 
shape of the model supports does not resemble the shape to be expected 
of actual supports. However, it is hoped that this model will provide 
insight into the scattering losses to be expected from an actual support 
structure. Since only order-of-magnitude estimates are expected from 
this model, we can use rather crude approximations of the mathe- 
matical expressions. The idea for using a naked fiber core held by a 
dielectric support structure originated with R. Kompfner. 


II. SCATTERING LOSS THEORY 


Our calculation of the scattering losses caused by the waveguide 
supports is based on the usual expansion of the field of the dielectric 
waveguide in terms of normal modes.! We restrict ourselves to TE 
modes of the slab waveguide whose only electric field component E, is 
tangential to the surface of the slab. The field is expanded in terms of 
normal guided and radiation modes of the slab waveguide 


By = 50.8, + | ao)8o)do. (1) 
fe 0 
The expansion coefficients in (1) are obtained from! 
L 
C, = o,f K pv(Z)e7* 89 Bw) eq , (2) 
0 


An analogous expression holds for g(p) with K,,, instead of A ,,. It can 
be shown [see eqs. (9.2-21) and (9.2-30) of Ref. 1] that the coupling 
coefficient K ,, is given by 

W EQ 2 * 

Ky =— (n(x) — nox) ]6,6,d2 . (3) 

4P J. 
no(x) 1s the refractive index distribution of the ideal waveguide, 
No(x) = ny, inside of the waveguide core and 29(X) = me outside of the 
core region. The index distribution (x) is defined as being n(a) = 11 
in the waveguide core and n(x) = n3 in the region of the waveguide 
support. ZL is the width of the support and P is the power carried by 
the incident mode. w and ¢ 9 are, respectively, the radian frequency of 
the field and the permittivity of vacuum. 
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The relative power loss from mode » to mode uw caused by the 
support is given by! 
L 
4 sin? (6, — — 
APy  |Cyl? (8, — By) 5 
= = | Ky]? ——_— (4) 
P; |C,| (8, _ By) 
In particular, this formula applies to the reflection loss R of the mode 


vy which is obtained by setting 8, = — 8,. The radiation loss of mode » 
is similarly obtained from the formula! 








L 
sin? (6, — B) — 


AP, rao) |2 ia 2 lal 
= ; Ig ic ge GG 
P, i ee a eee ae 


The integration variable in the integral on the extreme right-hand side 
was changed from p to 8. 

In order to evaluate the coupling coefficients, we need the field 
expressions for the guided and radiation modes only outside of the 
core region, because the integrand in (8) vanishes inside of the wave- 
guide core. We have for |x| >d 


Z 2yyP 
" U(n? = n2)8,(1 + yd)weo 








} 
eye 21a) | (6) 


Equation (6) for the guided TE modes follows from eqs. (8.38-12) and 
(8.3-18) of Ref. 1, with the help of (8.6-16). The parameters y, and 
x, are related to the free-space propagation constant k = 27/) in the 
following way 


ve = By — nak? (7) 
and 
Ps = nik? — B,. (8) 


The field of the TE radiation modes is obtained from eqs. (8.4-4), 
(8.4-9), and (8.4-18) of Ref. 1. 


p*k?P b 
B= jo 
2Q7we€o| 8 | (p? cos? cd + o? sin? od) 


Co 
. (co od — 4—-sin od) ete (lz|—4) 
p 


+ (cos od + ne sin cd) eo - (9) 
p 
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The coefficient for coupling between the guided modes » and u follows 
from (38) and (6) 


Na — Na (y7u)} 

= KyKy\ Vv ps 

Kyat sel gy 
mn = ns (Vy ++ VwL| BB u| (1 a yd) (1 lr ud) |} 





The power reflection coefficient FR, for the vth mode is obtained from 
(4) and (10) by setting 6, = — 8, «, =, andy, = vy 


(ns — no), sin? BL 1) 
‘ 4(n® — n2)*B4(1 + yd)? 
It is interesting to consider the reflection coefficient in the limit y, = 0 
at the cutoff point of the guided mode. We obtain from (11) with 
B, = nok and x2 = (nz — nb)k? 


(ng — nz)? 
R= a (4 sin? nokL)- (12) 
2 


The sine factor (multiplied by 4) on the right-hand side of this equation 
describes the interference between the reflection from the front and back 
surface of the support structure. If we omit this factor we obtain the 
reflection from only one of the two interfaces in the form 


poe (ng a ny)? _ (nz + ne)? (nz — 2)? 


7 16n' — (2n2)? (2nz)? U3) 


v 


Comparison of eq. (13) with the correct expression for the power 
reflection coefficient from a dielectric interface [see eq. (1.6-36) of 


Ref. 1 
etd N3 — No\? 
R = ( ) (14) 


N3 + Ne 





provides an indication of the accuracy of our approximation. Instead 
of solving the infinite system of coupled equations, we obtained the 
coefficient of the reflected wave and hence the power reflection coeffi- 
cient by considering only the incident and the reflected wave alone. 
The resulting equation (13) agrees with the correct equation (14) in 
the limit (n3 — 2) — 0. The approximate solution of the infinite equa- 
tion system is thus a good approximation only for small index differ- 
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ences. However, it is certainly quite satisfactory as an order-of- 
magnitude estimate. 

The coupling coefficient for mode » and a radiation mode charac- 
terized by the parameter p follow from (3), (6), and (9) 


2 2 
N3 — Ne 

a a ae } 
nem) 


kx,(y,)'p(y cos od — o sin od) 


ee 

(8° — 6*)[r|68,|(1 + 7,d)(p? cos? od + o? sin? ad) |} ve 
The radiative power loss caused by one support follows from (5) and 
(15). The integral cannot be solved exactly. Instead of resorting to 
numerical integration, I worked out an approximate solution which 
holds only as an order-of-magnitude estimate. 


AP — 2 (m3 — n2)* Ie» sas 
P 3a (n? — n2)# BL + ved)(B, + nak)*(8, — nok)? 


The radiation loss depends on the width L of the support structure. 
For large values of L/, the width dependence disappears. The approxi- 
mation (16) holds in this limit. Comparison of (16) with a few sample 
values of the numerical integration has shown that this approximation 
can depart from a few percent to as much as 50 percent from the value 
obtained by integrating (5). 


III. DISCUSSION AND NUMERICAL RESULTS 


The radiation losses caused by the waveguide supports are par- 
ticularly high for single-mode operation. Our discussion is thus directed 
towards multimode applications. However, Fig. 2 shows the depen- 
dence of the radiation loss on the width L of the supports for a single- 
mode case. It is apparent that the interference of radiation originating 
at the front and back surface of the support slab causes the radiation 
loss curve to oscillate. These oscillations die out with increasing width 
of the supports. The parameter V, that is used to label the curve in 
Fig. 1 and all remaining figures, is defined by 


V = (ny — nq) kd. (17) 


V is a parameter that combines frequency, slab width, and refractive 
index difference. Its values determine the number of modes that can 
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Fig. 2—Relative scattered power as a function of support thickness L. The dotted 
line represents eq. (16), which holds approximately for L >. 


propagate on the waveguide. The lowest-order even TE mode has no 
cutoff, the cutoff value of V for this mode is thus V, = 0. The cutoff of 
the next mode, the first odd mode, is given by V, = 7/2 = 1.57. In 
general, we obtain the cutoff condition for all the even and odd slab 


modes from the formula! 
fs 


Ve=v— 18 
an (18) 


For even modes, v assumes the values 0, 2, 4, etc.; for odd modes we 
have v = 1, 3, 5, ete. 

The solid curve of Fig. 2 was obtained by numerical integration of 
(5) and (15). The dotted line also shown in the figure results from the 
approximation (16). This approximation holds for L — © and is thus 
independent of the variable n2L/X. 

The following figures apply to multimode operation and show the 
radiation losses caused by one set of waveguide supports as a function 
of the mode angle 6. Hach guided mode can be decomposed into two 
plane waves, the propagation vectors of which form angles +6 and 
—6@ with the waveguide axis. The angle @ is defined by 


v 





cos 6, = (19) 


nik 
The mode angles assume discrete values 6, corresponding to the discrete 
values 8, of the propagation constant of the guided modes. Figures 3 
through 8 show 6 as a continuous variable. It is important to remember 
that only certain discrete values of 6 are allowed. The number of modes 
that exist below a certain value of 6 is approximately proportional to 
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Fig. 3—Relative scattered power (solid lines) as a function of mode angle 6. The 
waveguide core and the supports are In an air environment. The dotted lines represent 
the reflection loss of each mode. 


this angle in case of the slab waveguide but it is proportional to the 
square of 6 in case of the round optical fiber. For the slab, the number 
of TE and TM modes with angles smaller than a given value @ is 
approximately given by (@ < 1 is assumed) 


4 
No = — Onikd. (20) 
WT 


The total number of guided TE and TM modes that can be supported 
iS approximately given by 


4 
Nees = V. (21) 
Tv 
For the round fiber we have® (d = fiber radius) 
No = 3(On,4kd)? (22) 


modes with angles less than 6. The maximum number of modes that 
the fiber can support is® 
N max = $V?. (23) 
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Fig. 4—Scattering losses as a function of mode angle 6. The core index is assumed 
to be n, = 1.46, the support index is n3 = 1.3. Core and supports are immersed in 
an index matching liquid. The dotted line is the reflection loss for n3/nz = 1.01. 


Figure 3 shows the relative power loss as a function of mode angle 
6 for three different values of V. The refractive indices are chosen to 
represent a quartz fiber in air, 21/n2 = 1.46. It was assumed that the 
supports consist of a material with refractive index n3; = 1.3. The solid 
curves represent the radiation losses while the dotted curves show the 
reflection losses for the modes with mode angle @. It is apparent that 
the reflection losses are much smaller than the radiation losses. Mode 
conversion from each guided mode to all the other guided modes has 
not been considered. The maximum mode angle @ in this and all the 
following figures (with the exception of Fig. 4) coincides approximately 
with the right-hand edge of the graph. 

Figure 4 shows the radiation losses for a fixed value ni/n3 = 1.123 
of the ratio of core index to the index of the supports. This figure was 
drawn for the case that the core index is again n, = 1.46 and n3 = 1.3, 
but allows for the fact that an index matching liquid is used in an 
attempt to reduce the scattering losses. The ratio 3/72, that indicates 
the degree of index matching, is used as a curve parameter. The dotted 
curve is the reflection loss for n3/n2 = 1.01. The reflection losses for 
all the other index ratios are much smaller. The critical angle is 
different for each curve of this figure. 


214 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1973 










10-" np 1.005 Ls Vj 
0.141 ae we 
4 a 
wi 10-2 en | 
, Pa 
Lu 
QO 
; i ae 
~ lan 
; Rae 
107° 
0 0.02 0.04 0.06 0.08 0.1 012 0.14 
0 


Fig. 5—Relative scattering losses as a function of mode angle @, 1/n2 = 1.01, 
n3/Ne2 = 1.005. 
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Fig. 6—Same as Fig. 5, n3/n2 = 1.003. 
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Fig. 8—Same as Fig. 5, n3/n2 = 1.0001. 
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The remaining four figures are all similar to each other. All show 
the radiation losses as functions of mode angle for three different values 
of V and for different values of the ratio n3/ne, but fixed value 
3/N2 = 1.01. It is thus assumed that the index of the matching fluid 
remains the same, while the index of the supports is changed slightly. 
The reflection losses remain below the scale of all these figures. 

For an evaluation of the total radiation losses to be expected from 
a fiber held by many supports, it 1s necessary to know the total number 
of supports. Let us assume for simplicity that we have one support per 
centimeter. A fiber of 1-km length is thus held by a total of 10° sup- 
ports. In order to stay below a radiation loss of 10 dB/km, we must 
remain below the AP/P = 10-4 dB line of the figures. As a first crude 
approximation, we can assume that all modes, the loss of which 
remains below this line, are received at the end of the fiber while the 
modes that exceed this loss are lost. For a slab waveguide, the number 
of modes that can travel through the fiber is directly proportional to 
the mode angle 6, so that the ratio of transmitted to dissipated modes 
can be read off the horizontal axis. For the n3/n2 = 1.01 curve of Fig. 
4, modes with angles less than 0.26 radian stay below the 10-* line 
and can be considered as being transmitted through the fiber of 1-km 
length. This means that roughly half the modes are lost while half of 
them are being transmitted. However, if the waveguide is a round 
fiber (and if we accept the applicability of the loss curves for this case) 
we find that only one-quarter of all the modes can be transmitted while 
three-quarters of the modes and hence three-quarters of the power (for 
uniform initial power distribution) is lost. These estimates ignore mode 
conversion between the guided modes. The small values of the reflec- 
tion loss suggest that conversion between guided modes is reasonably 
small. However, some power is converted from modes with small angles 
to large angle modes so that the loss estimate made on the basis of 
the individual mode losses alone must be optimistic. 

A comparison of Fig. 4 and Fig. 5 shows clearly that it is advan- 
tageous to try to match the refractive index of the supports but keep 
the core index as large, compared to the surrounding liquid, as possible. 
The curves of Fig. 5 show much larger loss values because the difference 
in refractive index of core and index matching liquid is small, so that 
the fields extend farther into the liquid and thus are more effectively 
scattered by the supports. In order to achieve losses as low as those of 
the n3/n2 = 1.01 curve of Fig. 4, with a guide whose core-to-liquid 
index ratio is only 1.01, requires that the index ratio of support and 
matching liquid be better than 1.001. 
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IV. CONCLUSIONS 


The problem of radiative power loss caused by light scattering from 
the supports of an unclad fiber has been investigated. The study was 
based on a slab waveguide model with dielectric slabs used as supports. 
The study comes to the conclusion that tolerably low scattering losses 
are obtained only if the supports are made less visible to the wave by 
index matching with a suitable matching liquid. The matching liquid 
must itself have very low dielectric losses. It is important to make the 
index difference between core and supports as large as possible (with 
the core index being larger than that of the supports) and match the 
support index as closely as possible. For a core-to-support index ratio 
of ni/n3 = 1.123, an index ratio of n3/ne = 1.01 between supports 
and matching liquid is sufficient to allow one-quarter of all the modes of 
a round fiber to travel with acceptably low radiation losses. If the core 
index is more nearly equal to the refractive index of the supports, a 
much better index match for the supports is required. To achieve 
conditions comparable to the last example requires an index match of 
the supports of better than one tenth of a percent if n1/n2 = 1.01. It 
may be difficult to maintain such a good index match over the whole 
range of expected operating temperatures. 
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Passband Equalization of Differentially 
Phase-Modulated Data Signals 


By R. D. GITLIN, E. Y. HO, and J. E. MAZO 
(Manuscript received July 7, 1972) 


A new approach to the automatic equalization of differentially phase- 
modulated data signals passing through a linear distorting medium 1s 
presented. The equalizer, which 1s of the transversal filter type, operates 
in the frequency passband and contains two sets of taps—in-phase and 
quadrature. A tap-rotation property of the equalizer 1s used to establish 
an absolute phase reference at the equalizer output, and once this crucial 
step is taken, the passband equalizer output can be used to automatically 
(as well as adaptively) adjust the tap weights so as to minimize a mean- 
square distortion function. The resulting algorithm requires correlating an 
error signal with the tap voltages, and thus it is possible to use a structure 
similar to that employed when equalizing baseband PAM. 

The generality of the approach makes wt applicable to the equalization 
of any double-sideband data signal. 


I. INTRODUCTION 


Recently much attention has been given to high-speed synchronous 
data transmission via differential phase modulation and comparison 
detection. High-speed (above 2400 bits per second) transmission 
usually requires equalization! to compensate for the linear distortion 
introduced by the channel. This study will be concerned with the auto- 
matic and adaptive equalization of such data signals. At first glance, 
the nonlinear nature of the modulated signal would seem to preclude 
linear compensation. Upon closer examination it becomes apparent 
that a digitally phase-modulated signal, when resolved into in-phase 
and quadrature components, is linear in each component. Since this 
property is preserved after transmission through a linear medium, 
simultaneous linear equalization of each signal component becomes 
feasible. Due, however, to the purely quadratic nature of the channel- 
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dependent terms present in the differential detector output, a linear 
equalizer should’ precede the detector—hence the equalizer operates in 
the frequency passband. 

Resolving the channel output into in-phase and quadrature com- 
ponents suggests that satisfactory performance will be obtained if the 
equalized in-phase signal is approximately Nyquist and the equalized 
quadrature signal has most of its samples close to zero. This is accom- 
plished by choosing an appropriate cost function (of the equalized 
samples) and using an appropriately structured equalizer capable of 
minimizing this function. We will focus our attention on a mean-square 
cost function and on a synchronous tapped delay line (TDL) equalizer. 
The equalizer has two branches—an in-phase and quadrature branch-— 
which together perform the passband compensation. The delays are 
separated by a symbol interval, and each delay output is first multi- 
plied by a tap weight (each signal in the quadrature branch is also 
shifted by 90 degrees) and then added to the other delayed tap signals. 
The generality of the approach makes it applicable to the equalization 
of any double-sideband modulated data signal (e.g., combined ampli- 
tude and phase modulation, quadrature amplitude modulation). 

Though the differential detector output is independent of any 
absolute phase reference, it is convenient to select a cost function which 
depends on such a quantity. The awkwardness of this situation is re- 
solved by noting the tap-rotation property of the equalizer. Simply 
put, this means that any assumed reference phase angle can be 
“matched” by a rotation of the equalizer taps. This, in effect, permits 
the establishment of an arbitrary phase reference at the equalizer 
output. Equalization is accomplished with respect to this arbitrary 
reference while detection is done in an incoherent fashion. The tap 
adjustment algorithm is similar to that done in baseband PAM, te., 
tap signals are correlated with error signals. In fact, once the notion 
of a tap rotation is introduced, there is a convenient analogy to base- 
band PAM equalization. 

The basic system equations are indicated in Section II, and the 
equalizer cost function and structure are described in Section III. 
Section IV describes methods for adjusting the equalizer taps when 
an ideal data reference sequence is available and also when the adjust- 
ments are made using random data. In Section V we consider the 
effects of frequency offset and phase jitter on the equalizer. 


q In order that the mean-square distortion be a quadratic, and thus easily mini- 
mized, function of the in-phase and quadrature pulse samples. 
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Il. BASIC SYSTEM EQUATIONS 


In this section we indicate the response of a linear passband channel 
to a phase-modulated data signal. We also make some observations 
concerning the choice of phase reference that will be useful in our 
discussion of passband equalization. 

The channel input is a phase-modulated (PM) data signal. This 
signal is given by 


s(t) = > p(t — nT) cos (wet + On) (1a) 
= uin(t) COS wot — Vin(t) Sin wet (1b) 
= rin(t) cos (wet + Vin(t)), (1c) 


where 6, is the information symbol, 1/T is the symbol rate, w, is the 
carrier frequency, p(t) is the impulse response of the transmitter 
shaping filter,’ and 


Uin(t) = > p(t — nT) cos On (1d) 
Vin(t) = > p(t — nT) sin On (le) 


are respectively the in-phase and quadrature signal components. The 
envelope and phase of s;(¢) are given by 





ral) = Ve DFP] (af) 
Vi,(t) = tan7 1) (1g) 
Uin(E) 
By letting 
Gn = COS On (1h) 
b, = sin On, (11) 


we can write 


si(t) = [Xi anp(t — nT) ] cos we — [Xi bap(t — nT) ]sinwé. (1)) 


From eq. (1j), we can see that the bandlimited signal s;(é) is linear 
in both a, and b,.* As we have previously remarked, it is this linear 


' We assume the usual Nyquist shaping in the sense that p(k7') = 0, k ¥ 0. 

* Note that any double-sideband signal, such as quadrature amplitude modulation 
(QAM), can be written in the form of (1j). What we have to say in the sequel is 
sufficiently general to apply to any such double-sideband signals. 
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representation that makes linear compensation first plausible then 
possible. 

We note that the transmitter changes the value of @, once every T 
seconds, and that for differentially coded data, the customer infor- 
mation is the quantity 0, — 6,1. Applying s;(¢) to a bandpass channel 
with impulse response 


2F\(t) cos wt — 2F 2(t) sin wel, 
produces an output signal 


So(t) = u(t) cos wt — v(t) sin wel, (2a) 
where 
Uolt) = Fi) @uin(t) — F2(t) @vin(Z) 
volt) = Fo(t) @uin(t) + Fild) @vin (0), 


with @ denoting the convolution operation. To obtain a more compact 
representation, we define the in-phase channel pulse 


(2b) 


a(t) = pt) @Fi(), (3a) 
and the quadrature channel pulse 
y(t) = p(t) @F2(@). (3b) 


In terms of the data components 
Qn = COS 6,, 
b, = sin 6,, 
and the above notation, the channel output is given by 


So(t) = cos w,t[ >> anx(t — nT) — bray(t — nT) ] 


— sinw,f[ >) any(t — nT) + baxlt — nT)). (4) 


The above in-phase and quadrature representation clearly indicates 
the effect of the channel on the transmitted signal. A ‘‘good’”’ channel 
would be one which has a small quadrature pulse’ (ideally all its 
samples would be zero) and an in-phase pulse which is Nyquist. In 
any event, it is important to note that in passband the channel dis- 
tortion appears in a linear fashion. 


t If the channel has even amplitude symmetry and odd phase symmetry about the 
carrier, then y(t) will be identically zero for all ¢. 
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Since the carrier phase reference is arbitrary, the output signal can 
also be written as 


so(t) = cos (wt + O) [>> anxz(t — nT) — bray (t — nT) | 
— sin (wet + 6)[ D0 any(t — nT) + 6,245) (t — nT)] (5) 


where x‘ (t) and y‘*)(¢t) denote the in-phase and quadrature pulses 
with respect to the reference phase ¢. Comparing the right-hand sides 
(RHS) of (4) and (5) gives 


x?) (t) cos@ sin ¢ |[2(t) 6) 
eed ~ @ Cos La ( 

Equation (6) indicates that changing the phase reference by ¢ radians 
corresponds to rotating the in-phase and quadrature signal vector 
La(t), y(t) ] by this amount. We wish to emphasize the above in- 
terpretation, since we will want to consider the effect of such a rotation 
on several quantities of interest. We indicate the incoherent nature of 
the system by including an arbitrary phase angle when describing the 
channel output, i.e., the absolute phase of the received signal is un- 
known to the receiver. For incoherent reception, it is common to use 


differential coding of the data and differentially coherent (comparison) 
detection. We will assume this mode of transmission and detection. 


III. MEAN-SQUARE PASSBAND EQUALIZATION 


We wish to develop an adaptive equalizer which is capable of com- 
pensating for the linear (in-phase and quadrature) channel distortion, 
while keeping the equalizer adaption (settling) time small. ‘“Mean- 
square’ adaptive equalization! has been successfully used in systems 
using linear modulation to meet the above objectives. Our approach 
will be to extend, with appropriate modification, this technique to the 
problem at hand. After specifying the equalizer structure, we select 
an appropriate cost function and indicate how the equalizer processes 
the received random signal to adapt to the optimum configuration. 
We begin by considering the output of the differential detector. 


T We reserve the right to suppress the phase reference when there is no possibility 
of onfusion. Of course the phase reference can be chosen as convenience dictates. 
+ Clearly the same relation also holds for the sampled values of x(t) and y(t). 


224 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1973 


3.1 The Detector Output 


The equalizer is to be compatible with the standard comparison 
detector! which is shown in Fig. 1. It is straightforward to show that 
the detector outputs at the nth sampling instant are! 


I(nT + to) = 3» » cos (6; — 0;)[%n—-in—j-1 + Yn—-iY n—j-1] 


a Jj 


+ a 2 sin (0; —-0;)[ tances — YusFaj-1.) (7) 


v j 


QO(nT + to) = > ee COs (0; We 0;)[ Yn—iLn—j-1 a Ln—iYn—j-1 | 
tj 


+ 2d ps sin (6; =< 0;)[ ni n—j-1 + Yn—iY n—j_-1] (8) 


where 

a(n? +--4,) = 2, 
and 

y(nT + to) = Yn- 


We note that the intersymbol interference at the detector output is of 
a quadratic nature and thus is quite different than that encountered 
in linear modulation. It should also be noted that the channel-depen- 
dent terms in (7) are inner products, and hence are independent of the 
phase reference. The terms in (8) are components of cross products 
and hence also independent of phase. This is a manifestation of the 
differential nature of the receiver, i.e., only quantities which are in- 
variant with respect to a change in phase reference can affect the 
detector output. 

If we are fortunate enough to have a distortionless channel, 1.e., for 
all integers n, 


Ln = Sno (9a) 
an = A), (9b) 
then detection can be simply accomplished by noting that 
I(nT + to) = AnOn—-1 + Onbn—1 = cos(@n — On_1) (10a) 
Q(nT + to) = ba@n-1 — Gaba. = sin(@n — On-1). (10b) 


Unfortunately, (9) will not be satisfied for an arbitrary channel. As 
mentioned above, the function of the proposed equalizer is to linearly 


t We have, of course, neglected noise and have suppressed the phase reference. 
We also wish to emphasize that a desirable sampling epoch ¢, should be determined. 
As usual, this is a problem in its own right, but is not considered further in this paper. 
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Fig. 1—Comparison (differential) detector. 


operate on the in-phase and quadrature samples (the 2,’s and y,’s) 
to produce an overall system response the samples of which (denoted 
by g,’s and h,’s) approximately satisfy eq. (9). The quality of this 
approximation is measured by the value of an appropriate cost func- 
tion. We will have more to say about this later. 


3.2 Equalizer Structure 


Since we would like to equalize both the in-phase and quadrature 
channel samples, a structure which combines two signals in quadrature 
is suggested. The equalizer, which is shown in Fig. 2, consists of a 
tapped delay line (TDL) where each of the 2N + 1 delay outputs is 
fed into two branches. In each branch the delay outputs are multiplied 
by tap weights (c; and d;) and then summed. The equalizer output 
q(t) is the sum of the upper output qi(é) and the Hilbert transform of 
the lower output q2(t). The 90-degree phase shift at the output of the 
lower branch provides the quadrature signal. Thus the taps are to be 
adjusted to minimize an appropriate (cost) function of the equalized 
in-phase and quadrature samples. We remark that a linear equalizer 
placed in the passband can compensate for linear channel distortion, 
while a linear equalizer placed after the comparison detector cannot 
hope to compensate for the quadratic distortion present in the detector 
output. Thus the equalizer precedes the detector and operates in the 
frequency passband. We now describe the equalizer operation in detail. 
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Fig. 2—The passband equalizer. 


We can write the equalizer output, q(é), ast 
q(t) = cos atl >) ang(é — nT) — > DAE — nT) ] 


— sinwt >) an.htt — nT) + dO bag — nT)] (11) 


where the in-phase and quadrature (equalized) signals, g(t) and h(t), 
are given by? 


> crnz(t — nT) — > dny(t — nT) (12a) 


g(t) 


> dnx{t — nT) + > cny(t — nT). (12b) 


h(t) 


We emphasize that s;(¢) and q(t) are defined with respect to the same 


tT Assuming that w, and 7’ are chosen such that w.7' is an integer multiple of 27. 
If for some reason this is not convenient (perhaps a particular carrier frequency is 
desired), then w,7’ can be chosen to be any convenient angle. Some additional book- 
keeping is then required at the receiver. 

+ The time reference has been taken at the center tap of the equalizer. 


PASSBAND EQUALIZER 221 


phase reference, but that the phase reference has not been indicated 
in the above equations. For a sampling epoch ¢,, we denote the equal- 
ized samples by 


go = 9 GT +t.), hy =h(9T + 1,). (13) 


Thus with the reference phase taken to be ¢, the equalized samples 
are related to the channel samples and tap weights by 


N N 
(¢) (¢) (¢) 
Gy. i= PDR 7 oe 2 AnYj—n (14a) 
_N _N 
N N 
(¢) (¢) (¢) 
_N _N 


Again we reserve the right to suppress the representation phase in 
writing (14). We now make some observations which indicate the 
effects of a change in reference phase and of a “‘tap-rotation’”’ on the 
equalized samples. 

Denoting the channel phase reference by ¢, we can write the equalizer 
output as 


q(t) = cos (wt + O)[X ang(t — nT) — 2 baht — nT) I 
— sin (wet + O)L 2 anh (t — nT) + 2 brag (t — nT) J, (15) 


and it 1s easy to see that 


ae cos¢@ sin 4 ed 16) 
h‘(t) —sind@ cos¢ILh(t) 

Thus we again observe that changing the reference phase by ¢ radians 
corresponds to rotating the in-phase and quadrature signals (and 
samples) by this amount, i.e., each two-tuple (g@, h®) is rotated by 


@ radians. Suppose that while keeping the phase reference fixed, we 
transform the taps (c; and d;) to new values (¢; and d;) via the rotation 


C; cos @ sin 6][ ce; 
wel . i | f= = N60, 1D 
d; —sin @ cosé d; 


We can now express the equalized samples 2; and /;; (corresponding 
to ¢; and d;), using (14) and (17), as 


; cos @ sin @ ; | 
Hel | i i (Sas N sey. <8) 
h; —sin 6 cos 6JLh; 
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Thus rotating the equalizer taps by @ will rotate the equalized samples 
by the same amount.’ It is important to note that the original and 
(tap) rotated equalized samples are defined with respect to the same 
phase reference (¢). This is in contrast to the rotation introduced 
when the reference phase is changed. We will use these notions in 
describing the equalizer adjustment algorithm, in fact, the key to the 
simplicity of the equalization system lies in recognizing that the 
possibility of tap-rotation allows us to fix the reference phase and 
still keep a rotational degree of freedom between all (g;, h:) pairs. The 
reference phase will be suppressed in the sequel. 

It remains to choose an appropriate cost function (of the equalized 
samples) and to describe an iterative procedure for determining the 
tap weights that minimize this function. 


3.3 The Cost Function : 


The cost function should be such that after minimization, the 
equalized channel is approximately ideal (i.e., the intersymbol inter- 
ference has in some sense been suppressed). There are, of course, 
many such functions—thus our selection is influenced by additional 
considerations. 

We begin by noting that one possible task to set for the equalizer 
would be the minimization of 

D=Dgt Chi (19) 
1 40 1 0 
subject to appropriate constraint on go and ho.* We observe that since 
D is independent of an equalizer tap rotation, the minimum value of 
D is the same for any constraint that satisfies 


ea (20) 


including the choices 
Jo = 1, ho = 0. (21) 


Clearly a choice such as (21) fixes the reference phase used to describe 
the equalizer output, and a “‘rotation”’ of the equalizer taps can assure 
that it is met. Ultimately, the information needed to adjust the 
equalizer taps will indeed be obtained with respect to a particular 


' It should be clear that if there were no d, tap, we could not perform this rotation. 

* We note that this cost function would also be appropriate for a data signal which 
employs combined AM and PM. In fact, it is a simple matter to extend all of our 
results (i.e., the equalizer structure and algorithm) to this signaling format. 
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phase reference [in fact, we use the one implied by (21) ]. Hence the 
tap adjustments will be made using coherently obtained information. 

An additional interpretation of (21) is that the tap vectors are 
always required to lie on the hyperplanes go = 1, ho = 0. Unfortu- 
nately, minimizing a cost function subject to this (hard) constraint 
cannot be done in real-time. This is due to the awkward requirement 
that the tap vectors be adjusted while constrained to lie on hyper- 
planes. For example, suppose we were using a steepest descent? 
algorithm to adjust the taps, then we would need the projected (onto 
the constraint hyperplanes) gradient of D at every iteration.’ At 
present it is not known how to easily generate such a gradient from 
the available circuit voltages. 

We circumvent this difficulty by imposing a quadratic penalty (a 
soft constraint) when the variables do not satisfy the constraints.’ 
Thus the equalizer will try to minimize the unconstrained mean-square 
distortion 


= Cg t C+ — ge) +h. (22) 
t 40 t 0 
It should be noted that for small distortion, the optimum D [subject 
to (21) ] and & will be essentially identical.? 

A useful consequence of the notions of phase-reference rotation and 
tap rotation will now be demonstrated. Suppose the phase reference 
has been chosen and is fixed at this value. Then 6, which can be re- 
written as® 


6= D gith) +1— 2p (23) 
all z 

clearly depends upon the chosen phase reference, since go does. How- 
ever, since we have shown that a simple tap rotation can be used to 
effect any desired phase, we are guaranteed that the minimization of 
& will be over all reference phase angles. That is, while (23) is in general 
not reference-phase invariant, the minimum of (28) over all tap settings 
is reference-phase invariant. A further manifestation of the tap- 
rotation property is the following necessary condition for the optimum 
tap setting: The optimum tap setting is such that ho = 0. To see that 
this must be the case, suppose that the taps have settled down and 


t A more detailed discussion of this point is available in Ref. 4. 

+See Ref. 4 for a discussion of this type of subject for baseband PAM. 

§ Since the phase reference has been fixed, we do not indicate this quantity when 
writing the equalized samples. 
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Fig. 3—Illustration of the tap-rotation property. 


that ho ¥ 0. Since ©; 9? + h? is unchanged under a tap rotation, the 
equalizer taps could now be rotated so as to further minimize & by 
maximizing go (i.e., making hy = 0). The rotation of the vector (go, ho) 
into the vector (Vg + h2,0) is shown in Fig. 3; thus the required 
rotation is tan—!(ho/go). 

Recalling that a differential detector is only affected by phase- 
invariant quantities, we observe that if it were not possible to perform 
the above sort of minimization, then the utility of the structure and 
cost function would be suspect. 





IV. EQUALIZER TAP ADJUSTMENT 


Now that we have selected a cost function, we wish to demonstrate 
how the equalizer taps should be adjusted to minimize this function. 
We proceed in a manner reminiscent of mean-square equalization for 
baseband PAM-in fact it is shown that the tap increments are ob- 
tained by cross correlating error signals with delayed equalizer inputs. 


4.1 Deterministic Aspects 


We first wish to show that our criteria, 6, is an easily minimized 
function of the tap weights. We denote the tap vectors by 


cT = (c_y, «++, Coy +*°, €n) 


(24) 
dT = (d_y, +++, do, «+, dy) 


where the superscript denotes transpose. It is convenient to introduce 
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the channel correlation matricest 


Xe Dae (25a) 
EV dy = 2e Yn-iY n—j (25b) 
ie = 2 tas Saas) (25c) 


and the truncated channel vectors 


xt —_— (ty toe OG ++ *U_w) 
) ) (26) 
YS AYE © Yan Uy) 
We can now write the mean-square distortion, in matrix notation as 


rey ky 


= (7, a") | 
—~K X+YJ\d 


) = 207, -97 (|) 41. (27) 


A few comments are in order on the above representation. We first 
note that X and Y are the in-phase and quadrature channel correlation 
matrices. Thus if the quadrature component is zero, we need only the 
c taps. The cross-correlation between the in-phase and quadrature 
channel is measured by the skew-symmetric matrix K. By introducing 
the bandpass channel correlation matrix 


ast * + Y K | (28) 
—kK X+Y 

the augmented channel vector 
zt? = (27, —y?), (29) 

and the composite tap vector 
b? = (c?,d7), (30) 

the mean-square distortion can be expressed as 

6 = b?Ab — 2b7z2 + 1. (31) 


Conceptually at least, the optimum tap settings are given by setting 
the gradient of &, with respect to b, to zero. This gives? 


Ab —z=0 (32a) 


' LX ]:; denotes the 7jth element of the matrix X. | ae 
+ Positive definiteness of A is assumed. For the optimum tap setting it is easy to 
see that 6 = 1 — zTA7z. 
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or 
b* = Az (32b) 
as the desired tap settings. 

If a gradient (steepest-descent) algorithm is used to minimize’ & 
by iteratively adjusting b, then the speed of convergence of the 
algorithm (also referred to as the equalizer settling time) is determined 
by the matrix A. A simple bound on the speed of convergence is given 


in Ref. 4. 
1 Nin o 
(i) 
Amin Navex 


where Amax aNd Amin denote, respectively, the maximum and minimum 
eigenvalues of A, and &, is the mean-square error after the nth iter- 
ation. The ratio p = (Amax/min) Was computed for several voice-grade 
telephone channels (note that the closer p is to unity, the smaller the 
bound on settling time). Typically, the effect of the quadrature channel 
was small and the A matrix was diagonally dominant with p in the 
range from 2 to 8. Chang® has computed p for baseband PAM trans- 
mission, and comparing our numerical values with his indicates that 
for most channels the passband equalizer will settle rapidly. In other 
words, there are no special phenomena arising in the equalization of a 
bandpass channel which might lead one to expect a larger settling time 
than that observed for baseband channels. Though (32b) tells us what 
the optimum settings should be, we have not as yet indicated how 
these settings can be generated from available circuit voltages. 








Gn 
3 
IIA 


4.2 Real-Time Tap Adjustments Using an Ideal Reference 


Ultimately, the cost function actually minimized is determined by 
the ease with which the optimum tap settings can be obtained in real 
time. We first note that if, as assumed, A is positive definite, then & 
is a convex function of the tap vector b. Since, as is well known, a 
convex function has a single minimum, a gradient (steepest-descent) 
algorithm can be used to adjust the taps. If b“ denotes the tap vector 
after the kth adjustment, then the next tap setting is given by 


b&tD = bp — g®VE(b™), (33) 
where a“ is a gain (or step size) and V6& is the gradient of &. Note 


' The details of how one uses such an algorithm in this context are described in the 
next section. The reader not familiar with this type of algorithm should read Section 
4.2 before reading the present discussion. 
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that the algorithm ‘‘turns itself off’? only when the gradient is zero. 

Tor a careful choice of the sequence {a}, the above algorithm will 

converge to the optimum tap setting. Clearly the major problem 

associated with a real-time implementation of (83) is generating the 

gradient—or a good (statistical) approximation. To obtain an estimate 

of the gradient, we parallel the approach taken in baseband PAM. 
The gradient of 6 can be rewritten in component form, as 


ag 


oo oe se (Oitaas + hnYn=s) — Gi (34a) 
OC; n 

0&6 , 

ad. = (hnEn—i — GnYn—i) + Y-i Ls sd ++ +0, ea (34b) 


We remark that a reasonable statistical approximation (i.e., an esti- 
mate) of the true gradient would be one whose average value is the 
right-hand side of (34). To do this, we first use the sampled equal- 
izer output and an ideal data reference to generate a quantity, the 
average value of which is precisely our cost function. Once we have 
done this, we interchange the (linear) operations of ensemble averaging 
and differentiation to obtain the desired estimate. 
The equalizer output samples are given byt 


Q(jT + to) = qj = COS (weo + @)LX (Angj-n — Dnhj—n) ] 


— sin (woto + )[ > (Anhj—n -+ bngj—n) |, (35) 
and for perfect equalization these samples would be 


q; = 4; COS (welo + o) — 5; Sin (welo + b) = COS (wWebo +o; +). (86a) 


If we have an ideal reference (i.e., a; and b; are known at the receiver), 
then since the tap-rotation property permits the equalizer to obtain 
any phase, we have the liberty to choose ¢. In other words, we estab- 
lish a phase reference at the equalizer output, and we force (via the 
cost function penalty) the equalizer to find the phase angle of the 
vector (9%, h®) which minimizes the cost function. The equalizer 
does this by rotating the taps until h& = 0. For convenience we choose 
the reference phase, ¢, so that 


qj = (a; — b,), (36b) 


t Again, since the reference phase is arbitrary, we are free to choose ¢ as convenience 
dictates. 
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i.€., weto + @ has been taken to be 7/4. Note that we do not need an 
oscillator to generate the ideal reference. By analogy with baseband 
PAM, we next consider the mean-squared difference between the 
sampled equalizer output (q;) and the ideal (or desired) sequence 
(q'), i.e. EL(q; — 9). With 


a = Efa,] = E[b., (37a) 
the independence of the data sequences is used to show that 


EL(q — a) J = @ TX Gn thn) +1— 290] =a 8.  (87b) 


This is an extremely useful observation, since we now have available 
a quantity (q; — q3)°, the average value of which is proportional to the 
cost function. Interchanging the linear operations of expectation and 
differentiation allows us to write 


re] x 9 * 0q; 
ma EL(a; — 93) ] = 2k (a = qj) ae (38a) 

and ; 
| Gea) 128 fw — 4) at (38) 

ad, qj qj = Qj qj ad, 


Note that q; — q; is the error (at the equalizer output) at the jth 
sampling instant. It is easy to see that the differentiated terms can be 
written as 


0g; 
tue cos (wet + )[ >» AnXj—ii-n — a, On jin J 


Oc, n 
— sin (welo + O)[D0 AnYj-i_n + Do Dntj-i-n], (89a) 
and . | 


0g; 
= — cos (acto + o)[ > AnY jin + ys bata] 


— sin (weto + o)[ > Antj—i-n >»: OnYj—i—n |- (39b) 


We recognize that dq;/dc; is the jth channel output sample delayed by 
tT seconds, and thus is available at the zth delay element of the 
equalizer. We also note that 0q;/dd; is precisely 0q;/0c;, but with the 
carrier phase delayed by 7/2 radians, and thus is available at the out- 
put of the 7th delay followed by a 7/2 phase shifter. The estimated 
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gradient, which is given by the bracketed terms on the RHS of (88), 
can thus be obtained by correlating (muliplying and perhaps then 
averaging) the error signal (q; — q;) and the tap outputs.t As the 
equalization becomes better, the ideal reference can be replaced by 
decisions, i.e., the receiver operates in a decision-directed mode. This 
then gives a complete description of the passband equalization of the 
PM data signal when an ideal reference is available. 

At this juncture we wish to make the following remark. Suppose the 
channel is initially perfect [in the sense of eq. (9) ], but the phase 
reference is such that ho is not zero :* What does the equalizer do? We 
first observe that the equalizer will rotate the taps until ho Is zero, 
and thus lock in on the established reference phase. However, the 
differential detector output, which is insensitive to such a rotation, 
will have been error free from the outset. 


4.3 A Dectsion-Directed Tap Adjustment Procedure 


Suppose an ideal reference has been used to start-up the equalizer. 
It is desired that the equalizer then be capable of using decisions in 
place of an ideal reference. This mode of operation has been dubbed 
adaptive or decision-directed equalization.' The correction term needed 
for the tap adjustment algorithm is the product of the equalizer error 
signal (q; — q;) and a delayed equalizer input. For decision-directed 
operation, q; (which is a; — b;) would be replaced by the decided value 
of q;. For four-phase operation, a; and b; are binary, thus a; — b; 
= cos(¢; + 7/4) takes on the values —0.707 and 0.707. Thus g; can 
be obtained as the output of a threshold device (with a threshold at 
Zero). 


V. FREQUENCY OFFSET AND PHASE JITTER 


In this section we consider the effect of two common transmission 
impairments, frequency offset and phase jitter, on the operation of the 
equalizer. 


5.1 Frequency Offset 


Frequently, the transmission media perturbs the received carrier 
frequency w, so that it differs from the transmitted carrier frequency 


Fe embers that to obtain dq;/dd; the ith tap signal has been shifted by 7/2 
radians. 

* This could come about, for example, by a phase-hit impinging on an otherwise 
ideal channel. 
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27/T by A Hz.' As a result of this frequency offset, each tap signal has 
a different carrier phase. To obtain a more detailed understanding of 
this effect on the equalizer, we introduce the following rotated equalized 


samples 
n( A cosnA sin nA] gn 
bord i a a) 
h,(A) —sinnA cosnAILh, 


baer cosmA_ sin wl need (41) 
hin, n(A) ~ L=sin mA cosmA hi(A) J ’ 

where g, and h,, are the equalized samples in the absence of frequency 
offset. 


It is straightforward to show that the equalizer output, at the mth 
sampling instant, can be written in terms of the above quantities as 


qQimT) = Dy Gnxalmn\ A): = 2, OnsalnalA): (42) 


and 


The equalized samples [which are precisely the Zm,n(A)’s and the 
hm.n(A)’s] are thus a cascade of these two rotations. The rotation 
described by (40) is time-invariant and amounts to presenting the 
equalizer with a slightly different channel which is obtained by rotating 
the vector (gn, hn) by nA degrees. By observing that 


2 2 2 2 
mr gn(A) + An{A) = gn + hn (43) 
go(A) = go, 
it is clear that this effect can be neglected since the cost function is 
‘invariant under such a transformation. 

The rotation described by (41) 1s teme-varying and indicates that at 
the mth sampling instant, each channel pair (gn(A),hn(A)) is rotated 
by mA radians. Since each channel pair is rotated by the same amount, 
an equalizer tap-rotation can compensate for this rotation, the only 
requirement being that the equalizer settling time be much smaller 
than the period of rotation 1/A. Thus if the settling time of the equal- 
izer is small, the equalizer automatically tracks the frequency offset 
by rotating each tap pair (c,,d;) at the offset frequency A. 


5.2 Phase Jitter 


Phase jitter is an additive random component @(t) often present in 
the phase angle of the channel output, and is characterized as a low- 


' Typically A is less than 0.5 Hz. 
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pass process with energy up to 150 Hz. In this section we give only a 
suggestive description of the effect of phase jitter on the equalizer 
operation, since a detailed study appears to be difficult and could well 
be the subject of a separate investigation. If one were to assume the 
jitter to be constant along the length of the equalizer, 1.e., 


6(mT — iT) = 6(mT) i= —N,-::,0,°:-,N, 


then at the mth sampling instant, each equalized sample pair (gn, hn) 
is rotated by @(mT) radians. Such an assumption could only be valid 
for the extreme low-frequency components of the phase jitter, of 
course, and these would be successfully ‘‘tracked”’ by slow to-and-fro 
motions of the tap settings. The higher-frequency components could 
not be followed by the equalizer, but might be modeled by noise sources 
at each tap causing fluctuations of the individual settings about the 
optimum values. In such a model, it would be realistic to include cor- 
relations among these fluctuations at each spectral component and one 
would expect the correlations to decrease with increasing frequency. 


VI. SUMMARY AND CONCLUSIONS 


An approach to the mean-square equalization of a differentially 
phase-modulated data signal has been presented. An equalizer, of the 
transversal filter type, has been proposed which operates in the fre- 
quency passband and contains two sets of taps—-in-phase and quad- 
rature branches. By exploiting the tap-rotation property of the equal- 
izer, a phase reference is established at the output of the equalizer. 
Among the manifestations of the tap-rotation property are the ability 
to control the equalizer by using coherently obtained output samples 
and the ability to track small amounts of frequency offset and phase 
jitter. The equalized output is used to automatically (as well as 
adaptively) adjust the equalizer taps so as to minimize a mean- 
square distortion function. The required operations (correlating an 
error signal with tap voltages) are those performed when equalizing 
baseband PAM. Thus much of the existing knowledge concerning the 
technology of baseband PAM equalization can be applied to the 
equalization of phase-modulated data signals. 

We make two final remarks. The equalizer structure and tap adjust- 
ment algorithm can be applied, with minor modification, to the 
equalization of any double-sideband modulated data signals, e.g., 
combined amplitude and phase modulation. As is the case with all 
equalization systems, the precise dynamic behavior of the proposed 
equalizer can only be studied, at present, by experiment. 
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Selectively Faded Nondiversity and Space 


Diversity Narrowband Microwave 


Radio Channels 


By G. M. BABLER 
(Manuscript received August 14, 1972) 


The spectral characteristics of nondiversity and space diversity narrow- 
band radio channels subject to multipath fading were esteumated from a 
detailed sampling of channel loss variations as measured on two vertically 
separated antennas. The data base for this analysis was obtained during 
a 93-day experiment in which the amplitudes of a set of coherent tones 
spanning a band of 83.55 MHz and centered at 6084.2 MHz were con- 
tenuously monitored. The most significant observations were: 


(2) For the nondiversity channel, the frequency selectivity of the received 
transmission loss generally exceeded linear and quadratic components 
(in frequency) of amplitude distortion for fade depths greater than 
30 ab. 


(22) For the diversity channel constructed by switching between a pair of 
narrowband channels recewed on two vertically spaced antennas, the 
frequency selectivity of the transmission loss was significantly reduced. 


I. INTRODUCTION 


During periods of multipath propagation on line-of-sight radio links, 
deep and selective fading can severely corrupt the desired uniform 
transmission characteristics of a narrowband microwave channel re- 
ceived on a single antenna. Fortunately, because the multipath pro- 
cesses are sensitive functions of space,! the likelihood of simultaneously 
encountering both deep and frequency-selective fading on two verti- 
cally spaced antennas is very small. Earlier studies? were directed 
toward quantification of the spectral and temporal behavior of the 
amplitude distortion occurring in a nondiversity narrowband radio 
channel derived from a single receiving antenna. In the present paper, 
we describe a second experiment undertaken in 1971 to extend the 
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nondiversity measurements and to examine the amplitude distortions 
occurring in a diversity channel derived from the narrowband signals 
received on a pair of vertically spaced antennas. 

The data used in this analysis were taken between June 16 and 
September 16, 1971 (93 days), and included the amplitude measurc- 
ment of 62 uniformly spaced, coherent tones spanning 33.55 MHz at 
6 GHz transmitted over a 26.4-mile radio path and received on two 
antennas of vertical separation of 19 feet 3 inches. The tone fields re- 
ceived on the two antennas were sampled five times per second and 
the results recorded whenever significant amplitude activity was 
occurring anywhere in either channel. 

The paper is organized in the following order: (¢) description of the 
experiment, (zz) a discussion of the occurrence of fading throughout 
the 1971 fading season, and an identification of the data base for 
analysis, (777) characterizations of the nondiversity, and (2v) diversity 
amplitude distortions. 


II, SUMMARY 
A condensed listing of the findings follows: 


(z) The fade depth distributions for individual tones had the 
expected slopes? of a decade of probability of occurrence per 10 dB 
change in fade depth with the lower antenna undergoing somewhat 
less fading. 

(iz) The statistical distributions of linear and quadratic amplitude 
distortion of the nondiversity channel exhibited slopes of a decade of 
probability per 10 dB change in distortion. For a narrow channel 
bandwidth of 19.8 MHz, the linear and quadratic distortion exceeded 
18.5 and 7 dB, respectively, for 10-* of the observation time. 

(4it) For the nondiversity channel, the amplitude distortion gen- 
erally exceeded second order for fade depths greater than 30 dB and 
contained higher-order components in frequency selectivity. 

(iv) The statistical distributions of linear and quadratic amplitude 
distortion of the diversity channels derived by threshold or comparative 
switching exhibited slopes of a decade of probability of occurrence per 
5 dB change in distortion. For a bandwidth of 19.8 MHz, and for 10~° 
of the time, a —30-dB threshold diversity channel experienced linear 
and quadratic distortions exceeding 10 and 2.8 dB, respectively, and 
a comparative diversity channel experienced distortions of 8.6 and 
2.2 dB, respectively. 

(v) Yor the diversity channels, the amplitude distortion generally 
exceeded zeroth-order for fade depths greater than 10 dB. 
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Fig. 1—Experimental arrangement on the Atlanta to Palmetto, Georgia, radio link. 


III. EXPERIMENTAL DESCRIPTION 
3.1 The Experiment 


The 1971 measuring arrangement shown in Fig. 1 was essentially 
the same as the previous year which has been described elsewhere.’ 
To review briefly, a flat and constant in time narrowband field of 62 
coherent tones spaced 550 kHz apart and centered at a radio frequency 
of 6.0342 GHz was radiated from a standard horn reflector antenna 
at a microwave radio relay tower outside of Atlanta, Georgia. The 
signal, after propagating 26.4 miles along a line-of-sight path, was 
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Fig. 2—Subset of tones transmitted and recorded serially on the horn reflector 
and dish antennas at Palmetto. 
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intercepted by two antennas: a standard horn reflector with a 1.25- 
degree half-power beamwidth and a 6-foot dish with a 2-degree half- 
power beamwidth separated 19 feet 3 inches apart and mounted on a 
microwave tower at Palmetto, Georgia. Both images of the narrowband 
tone fields were translated to an intermediate frequency where a tone 
receiver selected tones for a power measurement. The subset of trans- 
mitted tones measured and recorded are shown in Fig. 2; the amplitude 
quantization was 1 dB, the time quantization was 0.2 second. A 
multiple input data acquisition system (MIDAS) supervised the 
measurement and recording of amplitude levels and time as well as 
controlled the recording rates according to the tone activity. 

The reference levels for both tone fields were determined by statis- 
tical studies of tone amplitudes during midday, nonfading periods. The 
rms variation in the reference levels was less than the amplitude 
quantization size of 1 dB. 


IV. OCCURRENCE OF FADING ACTIVITY 


The fading activity of the tone fields received on both antennas was 
continuously monitored from June 16 to September 16, 1971, and 
recorded for almost all of the 93 days (8.0352 * 10° seconds). The 
recorded data base, written by MIDAS and stored on 34 magnctic 
tapes, was condensed to include all periods with any fading in excess 
of 10 dB with respect to midday normal. This process compressed the 
time span of the data base to 1.28 X 10° seconds. 
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Fig. 3—Occurrence of fading for midchannel tone 24 received on the horn reflector 
antenna throughout the 1971 measurement period. 
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Fig. 4—Occurrence of fading for midchannel tone 24 received on the dish antenna 
throughout the 1971 measurement period. 


The data for a selected set of tones in both narrowband channels 
was processed to determine the total time during which any tone was 
faded below signal levels. An overview of the fading activity measured 
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Fig. 5—Fade depth distributions for tones 3, 24, and 60 received on the horn 
reflector antenna for the month of August. 
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Fig. 6—Fade depth distributions for tones 3, 24, and 60 received on the dish 
antenna for the month of August. 


during the fading season in Georgia and lumped into month categories 
of a single midchannel tone, tone 24, received on the horn and dish 
antennas is shown in Figs. 3 and 4, respectively. The occurrence of mid- 
channel fading for both antennas is far less in the June and September 
periods than for July or August. In addition, we observed more mid- 
channel fading for both antennas at and below 30 dB in July than in 
August, but the relative occurrence of fading at the 20, 30, and 40 dB 
levels and intermediate levels was more evenly balanced in August 
than in July. Because we desire to have a set of fading samples uni- 
formly spread across all fade depths with no particular levels dominat- 
ing the transmission loss statistics on either antenna, the month of 
August was selected as the working data base to be used in the char- 
acterization of the frequency-selective fading occurring in both 
channels. 

The fade depth distributions for tones 3, 24, and 60 received on the 
horn and dish antennas are given in Figs. 5 and 6, respectively. The 
ordinate is the fraction of 2.6784 10° seconds (number of seconds in 
August) that the tones were faded the amount indicated on the 
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Fig. 7—Unconditional linear and quadratic distortion distributions for the 19.8- 
MHz nondiversity (horn reflector antenna) channel. 


abscissa. The amplitudes, V, of the tones received on both antennas 
faded below signal levels, £, according to the expected fade depth 
distribution relation? 


PV sL)=al’ for Ls 0.1 (20 dB) — (1) 


where fade depth A = — 20 log L and where, for the August period, 
the proportionality constants a (functions of the radio paths’ electro- 
magnetic environment) were found to be a = 0.2 for the horn antenna 
and a = 0.1 for the dish antenna. The fact that, throughout the 
month, the lower dish antenna experienced less fading than the horn 
antenna reflects the spatial sensitivity of the fading processes and not 
the small differences in antenna beamwidths. Calculations based on an 
empirical formulation of the occurrence of multipath fading* indicate 
for the Atlanta-Palmetto 26.4-mile path at 6 GHz, a= 0.27. The 
conclusion is that, whereas the fading activity observed on both 
antennas was somewhat less than expected, the data base for the 
horn antenna is improved over the previous year.’ 

A second and important observation is that the fade depth distribu- 
tions of all three tones on both antennas are essentially the same, 
indicating a nonpreferential amount of fading in the channels which 
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was not the case in the samples of the 1970 study. This more uniform 
set of fade samples across the narrowband channels, as well as the 
somewhat greater number of events, result in an improved under- 
standing of frequency-selective fading in single and pairs of narrow- 
band radio channels. 


V. THE NONDIVERSITY CHANNEL SELECTIVITY CHARACTERIZATION 


The degree of frequency selectivity which occurred in the nondi- 
versity channel was characterized by statistical distributions of linear 
(AA =A (fs) —A (fi)) and quadratic (A?A/2= (A (f1) +A (fs))/2—-A (f2)) 
amplitude distortion constructed from three uniformly spaced samples 
of transmission loss in dB (A (f1), A (fe), A (f3)). The distributions for 
the linear and quadratic distortions occurring across a 19.8-MHz non- 
diversity channel (the signal from the horn antenna) are shown in 
Tig. 7. The abscissa is the amount of distortion in dB, and the ordinate 
is the fraction of 2.6784 X 10° seconds (August). The distributions 
exhibit slopes of a decade of probability of occurrence per 10 dB 
change in distortion with the linear distortion AA = 18.5 dB and the 
quadratic distortion A?A/2 2 7 dB for 10-5 of the time (about 27 
seconds). The roll-off of the curves below 10~° is a too-few-samples 
effect. For high-performance microwave radio systems the smallest 
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Fig. 8—Conditional linear distortion distributions for the 19.8-MHz nondiversity 
channel. 
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Fig. 9—Conditional quadratic distortion distributions for the 19.8-MHz non- 
diversity channel. 


fraction of time for single-link consideration is approximately 107°, or 
a few tens of seconds per month. For larger fractions of time the data 
show excellent point-to-point consistency. 

The linear and quadratic distortion distributions conditioned on fade 
depth of the midchannel tone [A(f:)] are shown in Figs. 8 and 9, 
respectively. The time base for each curve is indicated by the fade 
depth. Although the distortion curves have some point scatter, the 
curve-to-curve placement is improved over the 1970 sample.? In 
addition, we note that at each fade depth the channel experienced 
more linear and less quadratic distortion in 1971 as compared to the 
1970 sample. Because there were more deep fading events in 1971 
than 1970, and because there was more linear distortion at each fade 
depth in 1971 than 1970, we can conclude, as shown in Fig. 7 of this 
paper and Fig. 12 of Ref. 2, that there was more linear distortion in 
the channel. No similar arguments can be made for the quadratic dis- 
tortion because of conflicting occurrences of distortion and fade depth 
for both years, although we do note that the 1971 channel suffered 
less quadratic distortion. The rate of growth of the linear and quad- 
ratic distortion which occurred in the nondiversity channel with 
increasing bandwidth and fade depth was essentially the same as the 
1970 data (and is shown in Fig. 16). 
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Higher-order selective effects were again studied by accumulating 
the distributions of the maximum amplitude difference, called 
MAX|ERROR|, between the observed channel loss samples A™(f) 
and a three-term power series quadratic approximation constructed 
from the linear and quadratic distortion parameters. Thus 


MAX|ERROR| = MAX | AMD) 


—_——/ 


Measured Samples 
of Channel Loss 


- [aca + SG) + (Ga) ih 


Quadratic Power Series Approximation 


Fig. 10 shows the occurrence of MAX|ERROR| events for all fade 
depths (an unconditional distribution) and Fig. 11 shows the condi- 
tional distributions of MAX|ERROR| events for individual fade 
depths. The unconditional distribution of MAX|ERROR| for a band- 
width of 19.8 MHz is approximately the same as the 1970 data for a 
bandwidth of 20.835 MHz. The conditional distributions indicate 
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Fig. 10—Unconditional distribution of the MAX|ERROR| for the 19.8-MHz 
nondiversity channel. 
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Fig. 11—Conditional distributions of MAX|ERROR| for the 19.8-MHz non- 
diversity channel. 


rapid growth in the difference between the three-term power series 
approximation and the actual channel loss at and below the 30-dB 
fade level. For fade depths greater than 30 dB, the observed ampli- 
tude-frequency selectivity structure of the nondiversity narrowband 
radio channel’s transmission loss exceeded linear and quadratic com- 
ponents (in frequency) of amplitude distortion by 2 dB as was observed 
previously in the 1970 narrowband data. 


VI. THE DIVERSITY CHANNEL SELECTIVITY CHARACTERIZATION 


One form of protection against deep and selective fading in a narrow- 
band radio channel during periods of multipath propagation is to de- 
rive a diversity channel from the narrowband signals received on a pair 
of vertically spaced antennas. I*or the deeper fades, the reduction in 
number and duration of fading events increases significantly, since 
deep fades rarely occur simultaneously on two vertically spaced 
antennas.’ Because amplitude distortion (linear, quadratic, and higher 
orders of structure) is present for only the deeper fades, the likelihood 
of simultaneously encountering deep fading with large degrees of 
amplitude distortion on two antennas is indeed rare. 

To demonstrate the variety and nature of the measured transmission 
loss as simultaneously observed in the space diversity pair of narrow- 
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Fig. 12—Time sequential plots of the measured transmission losses simultaneously 
occurring in the narrowband channels received on the horn reflector and dish antennas. 


band channels during the 1971 period, and to motivate the diversity 
channel construction algorithms, two fading periods are presented in 
Figs. 12a and 12b. Displayed is the simultaneous time sequential state 
of the channel transmission loss of the tone fields as measured on the 
horn and dish antennas. An arbitrary 0-dB transmission loss reference 
has been employed for figure compactness; the time between scans 
was 0.2 second. 

In event (a), the selectivity develops at the upper edge of the narrow- 
band channel received on the horn antenna and then moves out of 
channel. The narrowband channel as received on the dish antenna was 
simultaneously undergoing about 5 dB of linear distortion super- 
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Fig. 13—Diversity channel construction schemes. 


imposed on a broader-band nonselective fade. In event (b), the horn 
antenna suffers only slight nonselective transmission loss, while the 
dish antenna experiences a rapidly moving selective fade. Most of 
the fading events observed were similar to events (a) or (b) ; generally, 
when there was deep fading with appreciable amplitude distortion in 
one channel, the alternate channel was undergoing only shallow fading. 

In Section 6.1 we describe how narrowband diversity channels were 
constructed from the space-diversity measurements of the 1971 
period and the characterization of the observed diversity narrowband 
channel amplitude distortions. 
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6.1 The Diversity Channel Construction 


The objective of diversity construction is to derive a new narrowband 
channel undergoing less fading and consequently less amplitude dis- 
tortion. Two straightforward idealized switching algorithms (schemes), 
threshold switching and comparative switching, have been examined 
and are presented schematically in Fig. 13. 

In threshold switching, a single narrowband threshold detector 
measures the power at the center of the narrowband tone field [say 
at A(fe) ] received on the presently connected antenna, and only if 
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Fig. 14—-Time sequential plots of the nondiversity and diversity channel trans- 
mission losses. 
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Fig. 15—Unconditional linear and quadratic distortion for the 19.8-MHz, —30-dB 
threshold diversity channel. 


this power falls below some fixed threshold level, system logic directs 
a radio frequency (RF) switch to select the alternate antenna. If the 
midchannel power on the secondary antenna is (or falls) below the 
threshold level, the main antenna is reconnected. The second switching 
scheme studied, comparative switching, employs two separate re- 
ceivers which simultaneously monitor the midchannel power of the 
narrow bands received on both antennas and direct the RF switch to 
connect in the antenna with greater midchannel power. Both forms of 
channel construction, threshold and comparative switching, result in 
a newly derived diversity channel undergoing less fading and con- 
sequently less amplitude distortion. 

For the selective events displayed in Fig. 12, both switching schemes 
would result in the selection of the less selectively faded antenna. The 
resulting diversity channel would therefore be the dish in event (a) 
and the horn in event (b). To serve as an additional instructive ex- 
ample, in Fig. 14 is shown the channel activity of a diversity channel 
constructed from the nondiversity channels which are simultaneously 
undergoing selective fading. The threshold diversity channel is the 
result of 10 switches between antennas beginning and ending the 
displayed period on the horn antenna. Because both antennas’ mid- 
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channel power levels are below the —30-dB threshold level, the cycle 
time between nondiversity channels is 0.2 second. For this particular 
set of events, contrary to those of Fig. 12, the threshold diversity 
channel’s amplitude distortion is not significantly improved over either 
nondiversity channel because of the simultaneously occurring selective 
fading in both antennas. The comparative diversity channel, the result 
of a single switch between antennas, experiences less amplitude dis- 
tortion than either nondiversity channel. The events as indicated in 
Tig. 12 are much more frequent than the highly selective and rare 
events as displayed in Fig. 14. 

By using the real-time measurements of transmission loss simul- 
taneously occurring in both narrowband channels and computer pro- 
grams to simulate the switching schemes, the lincar and quadratic 
amplitude distortion distributions for the diversity channels were 
accumulated both for unconditional and conditional fade depths 
[conditioned on the diversity channel’s tone, A(f2)] as well as for 
different bandwidths. These results follow. 


6.2 The Threshold Diversity Narrowband Channel 


The unconditional linear and quadratic distortion distributions for 
the 19.8-MHz-wide threshold diversity channel are shown in Fig. 15. 
A threshold level of —30 dB was chosen because it is at and below 
this level that the amplitude distortion of a nondiversity channel 
exceeds second-order selectivity and it is these higher-order distortion 
events we wish to avoid by the antenna selecting process. For a band- 
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Fig. 16—Growth of linear and quadratic distortion with bandwidth for the non- 
diversity channel and the —30-dB threshold diversity channel for a fixed 1075 
fraction of the time. 
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Fig. 17—Conditional linear distortion distributions for the 19.8-MHz, —30-dB 
threshold diversity channel. 


width of 19.8 MHz, the linear and quadratic distortions exceeded 10 and 
2.8 dB, respectively, for 10-° of the time. Contrasting Figs. 7 and 15, 
we find for 10~° fraction of the time that the —30-dB threshold switch- 
ing procedure has reduced the linear and quadratic amplitude dis- 
tortions by 46 and 60 percent, respectively, over the distortions 
occurring in the nondiversity channcl. The distributions of distortion 
for the —30-dB threshold channel exhibit slopes of a decade of prob- 
ability of occurrence per 5 dB change in distortion. Although the data 
points below 10~° are less reliable because of fewer samples, the roll-up 
at the tail of the linear distortion distribution is real and represents a 
few events spanning a few seconds during the 2.6784 X 10° seconds 
of measurement where severe fading below 30 dB was simultaneously 
occurring in both narrowband radio channels. The distributions for 
other bandwidths of the threshold channel were similar in shape to 
those presented in Fig. 15; the rate of growth in distortion with 
bandwidth, as indicated in Fig. 16, was about half the nondiversity 
channel’s rate. The —30-dB threshold diversity channel was the 
result of 103 switches between antennas. 

The conditional distributions of distortion are shown in Figs. 17 and 
18. The distributions show the expected reduction in linear distortion 
as compared to the nondiversity channel, Fig. 8, for fades not exceeding 
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Fig. 18—Conditional quadratic distortion distributions for the 19.8-MHz, —30-dB 
threshold diversity channel. 


30 dB. The roll-up of the distribution for 30-dB fades is the result of 
two effects: depletion of shallow structure (small linear distortion) 
events due to the threshold switching action and the over-abundance 
of a few high structure events at the 30-dB level. The roll-off of the 
35-dB fade distribution is the result of too few samples. The distribu- 
tions of the observed quadratic distortion at the various fade depths, 
Tig. 18, show a more uniform curve-to-curve behavior and indicate 
approximately the same amounts of distortion at each fade depth as 
the nondiversity channel. The 35-dB fade depth curve again has in- 
sufficient samples. 

These results show that a diversity channel selected according to the 
—30-dB threshold scheme outlined above undergoes significantly less 
overall amplitude distortion and presents to the radio sytem a more 
uniform and desirable transmission loss characteristic. Higher-order 
structure of the threshold diversity channel will be discussed in Section 
6.4 


6.3 The Comparative Diversity Narrowband Channel 


The unconditional linear and quadratic distortion distributions for 
the 19.8-MHz-wide comparative diversity channel are shown in Fig. 
19. Contrasting Figs. 15 and 19, we observe that the comparative di- 
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Fig. 19—Unconditional linear and quadratic distortion distributions for the 
19.8-MHz comparative diversity channel. 


versity channel experienced less linear and quadratic distortion than 
the threshold diversity channel. For a bandwidth of 19.8 MHz, the 
linear and quadratic distortions of the comparative diversity channel 
exceeded 8.7 and 2.2 dB, respectively, for 10-° of the time. Again 
we note that the distributions for the diversity channel exhibit slopes 
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Fig. 20—Growth of linear and quadratic distortion with bandwidth for the com- 
parative diversity channel for 10 of the time. 
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Fig. 21—Conditional linear distortion distributions for the 19.8-MHz comparative 
diversity channel. 


of a decade of probability of occurrence per 5 dB change in distortion. 
The scatter and roll-up at the tail for the linear distribution represents 
several seconds of fading events at which time both channels were 
undergoing significant linear amplitude distortion. The distributions 
for other bandwidths of the comparative diversity channel were similar 
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Fig. 22—Conditional quadratic distortion distributions for the 19.8-MHz com- 
parative diversity channel. 
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Fig. 23—Definition of deviations between the observed channel loss (solid line) 
and the zeroth degree, midchannel level approximation (dashed line). 


in shape to those presented in Fig. 19 and again the rate of growth in 
distortion with bandwidth was about half the nondiversity channel’s 
rate as indicated in Fig. 20. The comparative diversity channel was the 
result of 5873 switches between antennas. 

The conditional distributions of distortion are shown in Figs. 21 
and 22 and show approximately the same amounts of amplitude dis- 
tortion as the nondiversity channel for fade depths not exceeding 30 dB. 
The distributions for the 30-dB level have been included for complete- 
ness but are less reliable because of too few samples. 


6.4 Error Analysis of the Higher-Order Distortion 


As previously discussed in Section V, the selectivity structure within 
a narrowband channel increases with fade depth, and for the non- 
diversity channel, the selectivity structure exceeds second order for 
fades in excess of about 30 dB. For the diversity channels constructed 
by threshold and comparative switching, we found that those fading 
events with large linear and quadratic amplitude distortion were 
significantly reduced. Clearly, if a nondiversity channel does not 
exceed second order for fades not in excess of 30 dB, a diversity channel 
not undergoing fades in excess of 30 dB will also not exceed second 
order. Thus, rather than compute and construct the statistical distribu- 
tions for the maximum amplitude difference between the observed 
channel loss and a three-term power series quadratic approximation as 
we did for the nondiversity study, we choose to redefine the errors as 
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Fig. 24—Unconditional distributions of MAX |DEVIATION| for the nondiversity 
and diversity channels of bandwidth of 19.8 MHz. 


the deviations, D, between the measured samples of transmission loss, 


A™(f;), and a one-term power series approximation constructed from 
the midchannel loss A (fz). That is, 


Dp = Df) = -A™() — ACh). (2) 
ye ee Pe ee 
Measured Midchannel 
Samples of Loss 
Transmission 
Loss 


Figure 23 shows the observed selectivity, the zeroth-degree, midchannel 
level approximation, and two of the deviations. The maximum of the 
inband amplitude deviation, called MAX|DEVIATION|, was 
monitored for both forms of diversity channels and is now presented 
in Fig. 24. For completeness and for comparison, the nondiversity 
channel’s MAX |DEVIATION]| distribution is included. Both forms 
of diversity channels show significant reduction in maximum deviation 
as compared to the deviation experienced by the nondiversity channel. 
The marked reduction in deviations for the diversity channels was a 
result of less deep fading events in the diversity channels. 
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Fig. 25—Conditional distributions of MAX|DEVIATION| for the 19.8-MHz 
nondiversity channel. 


The amplitude deviation distributions conditioned on the mid- 
channel fade depth for the diversity channels were similar to those 
for the nondiversity channel which are presented in Fig. 25. For fade 
depths greater than 10 dB, the observed amplitude selectivity structure 
of the diversity (and nondiversity) narrowband channels exceeded 
zeroth order by 2 dB. This is in agreement with the observation that 
both forms of diversity switching may significantly protect the radio 
channel from deep selective fading events but not from the less faded 
channels which experience linear amplitude distortion. 
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A New Optical Fiber 
By P. Kaiser, EK. A. J. Marcatiui, and 8. E. MILLer 


(Manuscript received November 20, 1972) 


Currently there is strong interest in optical fibers for use as a trans- 
mission medium, analogous to the use of coaxial or wire pairs in the 
low-frequency region. Most work is devoted to a fiber structure consist- 
ing of a central glass core surrounded by a cylindrical glass cladding 
having a slightly lower index of refraction. This in turn requires that 
the chemical composition of the core glass differs from that of the clad- 
ding glass, leading to undesired effects at the core-cladding interface 
and perhaps limiting the minimum fiber losses achievable. 

The Nippon Sheet Glass Company and the Nippon Electric Com- 
pany together have developed a fiber (which they call SELFOC) with 
an index of refraction decreasing parabolically from the fiber axis to 
its outer boundary. This fiber requires a continuous variation in chem1- 
cal composition from the fiber axis outward, with attendant complica- 
tions in the fabrication process. A related guide requiring a film very 
thin compared to the wavelength has Just been reported.} 

The unique property of the new fiber is that a viable, handleable 
transmission medium is created by a structural form that uses only a 
single low-loss material. 

The conception was stimulated by the findings of P. Kaiser, et al., 
who fabricated unclad round fibers and measured their spectral losses 
in up to 32-meter unsupported lengths.? Recently he found total losses 
as low as 2.5 dB/km at wavelengths near 1.1 um using selected samples 
of low OH content fused silica. Similarly low losses have been measured 
at 1.06 um in bulk fused silica by T. C. Rich, et al.* It appeared attrac- 
tive to use material of this kind without the need for modifying the 
composition to alter the refractive index as is necessary with conven- 
tional core-cladding fibers or with graded-index fibers. 

Figure 1 shows section views of two possible forms of the single- 
material (SM) fiber. The usefully guided energy is concentrated pri- 
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Fig. 1—Cross section of an SM fiber with (a) rectangular and (b) cylindrical core. 


marily in the central enlargement, shown rectangular in Fig. la and 
round in Fig. 1b. In Fig. la, the central body has dimensions w and h 
for single-mode operation or dimensions A and B for multimode opera- 
tion. There is an exponentially decaying field extending outward from 
the central member in the slab of thickness ¢; with appropriate spacing 
between the central enlargement and outer cylinder the guided-wave 
field at the outside surface can be made negligibly small and the fiber 
can be handled exactly as can the conventional core-cladding type fiber. 
Slab modes are possible on the supporting structure, but these are 
strongly coupled to the outer shell and are readily lost to the surround- 
ing medium. Not shown in the figure is the possibility of adding an ab- 
sorbing coating on the outer surface for avoidance of crosstalk in a 
multifiber cable. 

The SM fiber structure can have a single propagating mode for any 
supporting slab thickness ¢ and for any shape of the central enlarge- 
ment, provided the size of the central enlargement is properly chosen. 
Practically, though, ¢ must be limited in order to keep the slab dimen- 
sion s, and consequently the overall size of the guide, reasonably small 
and still have the exponentially decaying slab field at the outer support- 
ing cylinder small enough. 

Analysis has been carried out for both single- and multimode SM 
fibers of several geometries. A few of the results are abstracted here. 
Yor the rectangular-guide case, Fig. la, and ¢ > X, there will be a single 
propagating mode provided 
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Note that wavelength does not appear in this expression, correct to 
first order. More exact analysis shows that for ¢ = 4.89 um andh = 7.0 
pm, the limiting width w for single-mode operation is 7.07 um at 
4 = 1.0um and 6.94 um at \ = 0.6 um. In these structures the slab 
field decays by 1/e in 2.80 and 2.75 um at » equal to 1.0 and 0.6 um, 
respectively. 

The wave propagation effects of the slab support can be represented 
by a uniform-index side support having the same height as the core and 
an equivalent index n. = n.(1 — A.), where n, is the index of the 
w-by-h core. Then, from the equality condition of eq. (1), it can be 


Ss Z 


Arbitrarily small values of A; can be achieved by making w [and accord- 
ing to (1), also h and ¢] appropriately large. 

For the multimode rectangular guide the number of guided modes 
may be shown to be 





cr AB 1 (3) 
7 2 4 (=) 
20 
where 
Deco 
v= ~ Vn? — n?, (4) 


and n is the index of the unshaded region outside the dotted region of 
Fig. la, and A and B are defined in the figure. Note that the number of 
modes, eq. (3), is (to first order) independent of wavelength—a unique 
property. For all modes the field decays exponentially along the slab 
as noted above; for the highest-order mode the field penetration is the 
largest and decays by 1/e in a length l, where 


v2t 7 \? 
= yi +(E) m 


For the multimode SM fiber the equivalent full-height support has an 
equivalent refractive index n.(1 ~— Am) where 


‘ =. vA \? 6) 
--<(;-) 


The value of A» can be used to calculate numerical aperture, the tol- 
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I'ig. 2—Photographs of an experimental (a) multimode SM fiber and (b) single- 
mode SM fiber (top), with magnified core region (bottom). 


erable radius of curvature, and modal dispersion. We give here only 
the numerical aperture, | 


ss d 
N.A. = n.V2A, = —: (7) 
2t 


SM fibers intended to approximate the geometries shown in Fig. 1 
were drawn in an oxygen-hydrogen torch from 6.5 mm i.d., 10 mm o.d. 
fused quartz tubes containing thin, polished plates and small-diameter 
rods supported in the center of the tubes. Plates of about 0.2 mm thick- 
ness, and core rods of approximately 0.2 mm and 1 mm diameters, 
resulted in single- and multimode fibers, respectively, whose cross sec- 
tions are shown in Figs. 2a and b. 

For 15-um-core multimode fibers, a slab thickness ¢ varying between 
3 and 4 wm resulted in numerical apertures ranging between 0.11 and 
0.08 (A = 0.6328 um), which agree excellently with the predicted 
values of 0.106 and 0.079, respectively [see eq. (7) |. The h/t ratio of 
the single-mode guide was about 6.5 um/4 um, or 1.625, with the width 
w amounting to 5 um. 

The spectral losses of the SM fibers were expected to closely approxi- 
mate those of the unclad fibers drawn from the same material. Whereas 
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this was true for the general shape of the spectral loss curve which was 
determined between 0.5 and 1.15 4m, the minimum losses were gen- 
erally higher. For a 300-m-long, Suprasil 2 multimode fiber they 
amounted to 39, 50, and 28 dB/km at 0.66, 0.80, and 1.06 um, respec- 
tively. Lowest losses of a single-mode fiber having a slightly different 
geometry than that shown in Fig. 2b were 55 dB/km at 1.06 um. We 
believe that residual contamination of the preform elements is the 
source of the excess losses. 

Total scattering losses in the order of 7.5 dB/km at 0.6328 um demon- 
strate that the approximately 30 modes (3) of the multimode fiber are 
well guided and do not lose power into the surrounding cladding to any 
significant degree. 

Other applications of the SM-fiber principle appear promising. 
Active fiber guides can be created by putting the active material in 
the central core or by putting it in a liquid surrounding the central 
member. Integrated optical circuits can utilize the same structure. For 
example, for a core and slab of index 1.472, slab thickness ¢ = 0.98 um, 
and a surround index 1 percent less than 1.472, we find the single-mode 
limit at h = 1.10 um and w = 6.75 um at »\ = 0.63828 um; the field 
decays transversely in the slab by 1/e in 2.67 um. Thicker slabs allow 
larger w and h with single-mode guidance. In early research on optical 
integrated circuits, J. E. Goell observed wave propagation in a curved 
guide of the above general form, now understood as another verification 
of the principle of the SM fibers. 

The assistance of H. W. Astle in the fabrication of the SM fibers is 
gratefully acknowledged. 
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